A DEODHAR DECOMPOSITION OF THE GRASSMANNIAN AND THE 
REGULARITY OF KP SOLITONS 



(N 



YUJI KODAMA AND LAUREN WILLIAMS 



Abstract. Given a point A in the real Grassmannian, it is well-known that one can construct a 
£N| soliton solution ua(x, y, t) to the KP equation. The contour plot of such a solution provides a tropical 

, | approximation to the solution when the variables x, y, and t arc considered on a large scale and the 

Q time t is fixed. In this paper we use several decompositions of the Grassmannian in order to gain an 

understanding of the contour plots of the corresponding soliton solutions. First we use the positroid 
stratification of the real Grassmannian in order to characterize the unbounded line-solitons in the 
contour plots at y 2> and 1/ « 0. Next we introduce a refinement of the positroid stratification - 
the Deodhar decomposition of the Grassmannian — which is defined to be the projection of Dcodhar's 
decomposition of the complete flag variety. We index the components of the Deodhar decomposition 
of the Grassmannian by certain tableaux which we call Go-diagrams, and then use these Go-diagrams 
to characterize the contour plots of solitons solutions when t«0. Finally wc use these results to 
r \ show that a soliton solution u^(x, y, t) is regular for all times t if and only if A comes from the totally 

non-negative part of the Grassmannian. 
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1. Introduction 

The KP equation is a two-dimensional nonlinear dispersive wave equation which was proposed by 
Kadomtsev and Peviashvili in 1970 to study the stability problem of the soliton solution of the Korteweg- 
de Vries (KdV) equation [13]. The KP equation can also be used to describe shallow water waves, and 
in particular, the equation provides an excellent model for the resonant interaction of those waves. The 
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equation has a rich mathematical structure, and is now considered to be the prototype of an intcgrablc 
nonlinear dispersive wave equation with two spatial dimensions (see for example [25, 1, 9, 24, 12]). 

One of the main breakthroughs in the KP theory was given by Sato [30] , who realized that solutions 
of the KP equation could be written in terms of points on an infinite-dimensional Grassmannian. The 
present paper deals with a real, finite-dimensional version of the Sato theory; in particular, we are 
interested in solutions that are localized along certain rays in the xy plane called line-solitons. Such 
a soliton solution can be constructed from a point A of the real Grassmannian. More specifically one 
can apply the Wronskian form [30, 31, 11, 12] to A to produce a r-function TA{x,y,t) which is a sum 
of exponentials, and from the r-function one can construct a solution UA{x,y,t) to the KP equation. 

Recently several authors have studied the soliton solutions ua{x, y, t) which come from points A of the 
totally non-negative part of the Grassmannian (Grk. n )>a, that is, those points of the real Grassmannian 
Gr k n whose Pliicker coordinates are all non-negative [3, 17, 2, 5, 7, 19, 20]. These solutions are regular, 
and include a large variety of soliton solutions which were previously overlooked by those using the 
Hirota method of a perturbation expansion [12]. 

One of the main goals of this paper is to understand the soliton solutions UA(x,y,t) coming from 
arbitrary points A of the real Grassmannian, not just the totally non-negative part. In general such 
solutions are no longer regular - they may have singularities along rays in the xy plane - but it is 
possible, nevertheless, to understand a great deal about the asymptotics of such solutions. 

Towards this end, we use two related decompositions of the real Grassmannian. The first decomposi- 
tion is Postnikov's positroid stratification of the Grassmannian [27], whose strata are indexed by various 
combinatorial objects including decorated permutations and J-diagrams. Note that the intersection of 
each positroid stratum with (Grk, n )>o is a cell (homeomorphic to an open ball); when one intersects 
the positroid stratification of the Grassmannian with the totally non- negative part, one obtains a cell 
decomposition of (Grk. n )>o [27]. 

The second decomposition is a refinement of the positroid stratification which we call the Deodhar 
decomposition of the Grassmannian. This is obtained by projecting the Deodhar decomposition of the 
complete flag variety to the Grassmannian. Its components have explicit parameterizations due to 
Marsh and Rietsch [23], and are indexed by distinguished subexpressions of reduced words in the Weyl 
group. The components may also be indexed by certain tableaux filled with black and white stones which 
we call Go-diagrams, and which provide a generalization of J-diagrams. Note that almost all Deodhar 
components have an empty intersection with the totally non-negative part of the Grassmannian. More 
specifically, each positroid stratum is a union of Deodhar components, precisely one of which has a 
non-empty intersection with (Grfc,„)>o- 

By using the positroid stratification of the Grassmannian, we characterize the unbounded line-solitons 
of KP soliton solutions coming from arbitrary points of the real Grassmannian. More specifically, given 
A e Grk.n, we show that the unbounded line-solitons of the solution UA(x,y,t) at y <C and y >• 
depend only on which positroid stratum A belongs to, and that one can use the corresponding decorated 
permutation to read off the unbounded line-solitons. This extends work of [2, 5, 7, 19, 20] from the 
setting of the non-negative part of the Grassmannian to the entire real Grassmannian. 

By using the Deodhar decomposition of the Grassmannian, we give an explicit description of the 
contour plots of soliton solutions when f<0. More specifically, we provide an algorithm for constructing 
the contour plot of UA{x,y,t), which uses the Go-diagram indexing the Deodhar component of A. We 
also show that when the Go-diagram D is a J-diagram, then the corresponding contour plot at t <C 
gives rise to a positivity test for the Deodhar stratum stratum Sd- 

Finally we use our previous results to address the regularity problem for KP solitons. We prove that 
a soliton solution UA(x,y, t) coming from a point A of the real Grassmannian is regular for all times t 
if and only if A is a point of the totally non-negative part of the Grassmannian. 
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The structure of this paper is as follows. In Section 2 we provide background on the Grassmannian 
and some of its decompositions, including the positroid stratification. In Section 3 we describe the 
Deodhar decomposition of the complete flag variety and its projection to the Grassmannian, while in 
Section 4 we explain how to index Deodhar components in the Grassmannian by Go-diagrams. In 
Section 5 we provide explicit formulas for certain Plucker coordinates of points in Deodhar components, 
and use these formulas to provide positivity tests for points in the real Grassmannian. Subsequent 
sections provide applications of the previous results to soliton solutions of the KP equation. In Section 
6 we give background on how to produce a soliton solution to the KP equation from a point of the real 
Grassmannian. In Section 7 we define the contour plot associated to a soliton solution at a fixed time t, 
then in Section 8 we use the positroid stratification to describe the unbounded line-solitons in contour 
plots of soliton solutions at y » and yCO. In Section 9 we define the more combinatorial notions of 
soliton graph and generalized plabic graph. In Section 10 we use the Deodhar decomposition to describe 
contour plots of soliton solutions for t <C 0, and in Section 11 we provide some technical results on 
X-crossings in contour plots and corresponding relations among Plucker coordinates. Finally we use 
the results of the previous sections to address the regularity problem for soliton solutions in Section 12. 

2. Background on the Grassmannian and its decompositions 

The real Grassmannian Gr kn is the space of all fc-dimensional subspaces of R™. An element of Gr kiU 
can be viewed as a full- rank k x n matrix modulo left multiplication by nonsingular k x k matrices. In 
other words, two k x n matrices represent the same point in Gr k , n if and only if they can be obtained 
from each other by row operations. Let be the set of all fc-element subsets of [n] := {1, . . . , n}. 
For / G ('fe')j l e t Ai(j4) be the Plucker coordinate, that is, the maximal minor of the k x n matrix A 
located in the column set /. The map A ^ (Aj(A)), where I ranges over (^'), induces the Plucker 

embedding Gr k , n ^ WOO -1 . 

We now describe several useful decompositions of the Grassmannian: the matroid stratification, the 
Schubert decomposition, the positroid stratification, and a Deodhar-type decomposition. 

2.1. The matroid stratification of Gr k 

Definition 2.1. A matroid of rank k on the set [n] is a nonempty collection M C (^') of fc-element 

subsets in [n], called bases of M, that satisfies the exchange axiom: 

For any I,JeM and i e I there exists j € J such that (/ \ {i}) U {j} £ M. 

Definition 2.2. A loop of a matroid on the set [n] is an element i € [n] which is in every basis. A coloop 
is an element i e [n] which is not in any basis. 

Given an clement A e Grk, n , there is an associated matroid Ma whose bases are the fc-subsets 
I C [n] such that A/ (A) ^ 0. 

Definition 2.3. Let M C (^') be a matroid. The matroid stratum Sm is defined to be 

S M ={AE Gr k . n | A/(A) ^ if and only if I e M}. 

This gives a stratification of Gr k n called the matroid stratification, or Gelfand-Serganova stratification. 
The matroids M. with nonempty strata Sm ar e called realizable over K. 

2.2. The Schubert decomposition of Gr k ^ n . We now turn to the Schubert decomposition of the 
Grassmannian. First recall that the partitions A C (n — k) k are in bijection with /c-element subset 
/ C [n]. The boundary of the Young diagram of such a partition A forms a lattice path from the 
upper-right corner to the lower- left corner of the rectangle (n — k) k . Let us label the n steps in this 
path by the numbers 1, . . . ,n, and define / = /(A) as the set of labels on the k vertical steps in the 
path. Conversely, we let A(J) denote the partition corresponding to the subset /. 
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Definition 2.4. For each partition A C (n — k) k , one can define the Schubert cell fl\ by 

fix = {A € Grfe,„ | /(A) is the lexicographically minimal base of Ma}- 

As A ranges over the partitions contained in (n — k) k , this gives the Schubert decomposition of the 
Grassmannian Gr kn , i.e. 

Gr fej „ = ;J fl x - 

\<Z(n-k) k 

We now define the shifted linear order <j (for i e [n]) to be the total order on [n] defined by 
i <i i + 1 <i i + 2 <i • • • <i n <i 1 <i • • • <i i — 1. 
One can then define cyclically shifted Schubert cells as follows. 

Definition 2.5. For each partition A C (n — k) k and ie [n], we define the cyclically shifted Schubert cell 
fl\ by 

0\ = {A e Gr^„ | /(A) is the lexicographically minimal base of Ma with respect to <j}. 
Note that fl x = fl\. 

2.3. The positroid stratification of Gr k ^ n . The positroid stratification of the real Grassmannian 
Grk.n is obtained by taking the simultaneous refinement of the n Schubert decompositions with respect 
to the n shifted linear orders <j. This stratification was first considered by Postnikov [27], who showed 
that the strata are conveniently described in terms of Grassmann necklaces, as well as decorated per- 
mutations and 1-diagrams. Postnikov coined the terminology positroid because the intersection of the 
positroid stratification with the totally non-negative part of the Grassmannian (Grfc jTl )>o gives a cell 
decomposition of (Grk, n )>o (whose cells are called positroid cells). 

Definition 2.6. [27, Definition 16.1] A Grassmann necklace is a sequence X = . . . ,I n ) of subsets 
I r C [n] such that, for i e [n], if i £ Ii then I i+ i = (7, \ {i}) U {j}, for some j e [n]; and if i ^ I r then 
I i+ i = Ii. (Here indices i are taken modulo n.) In particular, we have |7i| = • • • = |7„|, which is equal 
to some k € [n]. We then say that X is a Grassmann necklace of type (fc,n). 

Example 2.7. 1 = (1257,2357,3457,4567,5678,6789,1789,1289,1259) is an example of a Grassmann 
necklace of type (4,9). 

Lemma 2.8. [27, Lemma 16.3] Given A e Gr fei „, let I{A) = (7 1; . . . ,7„) be the sequence of subsets in 
[n] such that, for i€ [n], Ii is the lexicographically minimal subset of (^') with respect to the shifted 
linear order <j such that Aj. (^4) ^ 0. Then 1(A) is a Grassmann necklace of type (k,n). 

If A is in the matroid stratum Sm, we also use Im to denote the sequence (7i, . . . , 7„) defined above. 
This leads to the following description of the positroid stratification of Gr k ^ n . 

Definition 2.9. Let I = (I\, . . . , 7„) be a Grassmann necklace of type (k,n). The positroid stratum Sx 
is defined to be 

Si = {A g Gr fc ,„ | 1(A) = I}. 

Remark 2.10. By comparing Definition 2.9 to Definition 2.5, we see that given a Grassmann necklace 
I = (7i, . . . ,7„), 

n 

Si = Pi ^x(h) ■ 
i=i 

In other words, each positroid stratum is an intersection of n cyclically shifted Schubert cells. 
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Definition 2.11. [27, Definition 13.3] A decorated permutation 7T = (n,col) is a permutation tt G S n 
together with a coloring function col from the set of fixed points {i \ n(i) — i} to {1, —1}. So a decorated 
permutation is a permutation with fixed points colored in one of two colors. A weak excedance of w is a 
pair (i, 7Tj) such that either 7r(i) > i or 7r(i) = i and col(i) — 1. We call i the weafc excedance position. 
If 7r : (i) > i (respectively 7r(i) < i) then (i,7Ti) is called an excedance (respectively, nonexcedance). 

Example 2.12. The decorated permutation (written in one-line notation) (6,7,1,2,8,3,9,4,5) has no 
fixed points, and four weak excedances, in positions 1, 2, 5 and 7. 

Definition 2.13. [27, Definition 6.1] Fix fc, n. If A is a partition, let Y\ denote its Young diagram. A 
J-diagram (A, D)k^ n of type (k, n) is a partition A contained in a k x (n — fc) rectangle together with a 
filling D : Y\ — > {0, +} which has the J-property: there is no which has a + above it and a + to its 
left. 1 (Here, "above" means above and in the same column, and "to its left" means to the left and in 
the same row.) 

In Figure 1 we give an example of a J-diagram. 
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k = 6, n=16 

K = (10, 9, 9, 8, 5, 2) 



FIGURE 1. A Le-diagram L = (A,D) fe ,„. 

We now review some of the bijections among these objects. 

Definition 2.14. [27, Section 16] Given a Grassmann necklace I, define a decorated permutation w = 
7T (X) by requiring that 

(1) if = (h \ {»}) U {j}, for 3 ? i, then tt : (j) = i. 2 

(2) if 2j_|_i = Ii and i G Zj then 7r(i) = i is colored with col(i) = 1. 

(3) if = and i ^ Ii then 7r(z) = i is colored with col(i) = — 1. 
As before, indices are taken modulo n. 

If 7T = 7r(Z), then we also use the notation S v -- to refer to the positroid stratum St- 

Example 2.15. Definition 2.14 carries the Grassmann necklace of Example 2.7 to the decorated permu- 
tation of Example 2.12. 

Lemma 2.16. [27, Lemma 16.2] The map T — > 7r(Z) is a bijection from Grassmann necklaces I = 
(Ji, . . . , I n ) of size n to decorated permutations 7r (I) of size n. Under this bijection, the weak excedances 
of 7T (I) are in positions I\ . 

Remark 2.17. Use the notation of Lemma 2.16. It follows from the definition of the positroid stratifica- 
tion that if A G Si is written in row-echelon form, then the pivots are located in position I\ . It follows 
from Lemma 2.16 that the pivot positions coincide with the weak excedance positions of 7r(I). 



This forbidden pattern is in the shape of a backwards L, and hence is denoted J and pronounced "Le." 
2 Actually Postnikov's convention was to set TT'(i) = j above, so the decorated permutation we are associating is the 
inverse one to his. 
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2.4. Irreducible elements of GVfc T 



Definition 2.18. We say that a full rank k x n matrix is irreducible if, after passing to its reduced row 
echelon form A, the matrix A has the following properties: 

(1) Each column of A contains at least one nonzero element. 

(2) Each row of A contains at least one nonzero element in addition to the pivot. 

An irreducible Grassmann necklace of type (k, n) is a sequence I — . . . , /„) of subsets I r of [n] of 
size k such that, for ie [n], = (h \ {i}) U {j} for some j ^ i. (Here indices i are taken modulo n.) 
A derangement n — (jti, . . . , 7r„) is a permutation n € S n which has no fixed points. 

In the language of matroids, an element A e Sm is irreducible if and only if the matroid M. has 
no loops or coloops. It is easy to see that if A is irreducible, then 1(A) is an irreducible Grassmann 
necklace and 7r(Z) is a derangement. 

3. Projecting the Deodhar decomposition of G/B to the Grassmannian 

In this section we review Deodhar's decomposition of the flag variety G/B [8], and the parameteriza- 
tions of the components due to Marsh and Rietsch [23] . We then give a Deodhar decomposition of the 
Grassmannian by projecting the usual Deodhar decomposition of the flag variety to the Grassmannian. 

3.1. The flag variety. The following definitions can be made for any split, connected, simply con- 
nected, scmisimple algebraic group over a field K. However this paper will be concerned with G = 
SL„ = SL„(IR). 

We fix a maximal torus T, and opposite Borel subgroups B + and B~ , which consist of the diagonal, 
upper-triangular, and lower-triangular matrices, respectively. We let U + and U~ be the unipotcnt 
radicals of B + and B~ ; these are the subgroups of upper-triangular and lower-triangular matrices with 
l's on the diagonals. For each 1 < i < n - 1 we have a homomorphism : SL 2 — > SL„ such that 

A \ 



a b 

c d 



a b 
c d 



e SL ri 



V V 

that is, 4>i replaces a 2 x 2 block of the identity matrix with ^ ^ . Here a is at the (i + l)st 

diagonal entry counting from the southeast corner. 3 We use this to construct 1-parameter subgroups 
in G (landing in U + and U~ , respectively) defined by 



Xi(m) = 4> l 



1 m 
1 



and yi(m) = cf) l 



1 
m 1 



where m G 



The datum (T, B + , B , Xi, yi \ i € /) for G is called a pinning. 

Let W denote the Weyl group Nq(T)/T, where Nq(T) is the normalizer of T. The simple reflections 

Si e W are given explicitly by Sj := s-iT where s\ := fa ( ^ r\ ) anc ^ an y w can be expressed 



V 1 

as a product w = s il s i2 . . . s im with m = £(w) factors. We set w 
W = & n , the symmetric group on n letters, and is the transposition exchanging i and i + 1. 



SjjS^ . . . s im . For G — SL„, we have 



3 Our numbering differs from that in [23] in that the rows of our matrices in SL n arc numbered from the bottom. 
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We can identify the flag variety G/B with the variety B of Borel subgroups, via 

gB^g-B+ :=gB+g-\ 
We have two opposite Bruhat decompositions of B: 

B= \_\ B+w-B+= \_\ B-V-B+. 

Note that B~v ■ B + = ]R £ ("' ) _£ ( 1 '). The closure relations for these opposite Bruhat cells are given by 
B~i>' ■ B + C B~v ■ B+ if and only if v < v' . We define 

Tl VyW := B+w ■ B+ n B~v ■ B+, 

the intersection of opposite Bruhat cells. This intersection is empty unless v < w, in which case it is 
smooth of dimension £(w) — £(v), see [15, 22]. The strata 1Z V>W are often called Richardson varieties. 

3.2. Distinguished expressions. We now provide background on distinguished and positive distin- 
guished subexpressions, as in [8] and [23]. We will assume that the reader is familiar with the (strong) 
Bruhat order < on the Weyl group W — 6„, and the basics of reduced expressions, as in [4]. 

Let w := Sj x . . . s im be a reduced expression for w G W. We define a subexpression v of w to be a word 
obtained from the reduced expression w by replacing some of the factors with 1. For example, consider 
a reduced expression in 64, say S3S2S1S3S2S3. Then s 3 s 2 1 s 3 S2 1 is a subexpression of S3S2S1S3S2S3. 
Given a subexpression v, we set V(k) to be the product of the leftmost k factors of v, if k > 1, and 
f(o) = 1. The following definition was given in [23] and was implicit in [8]. 

Definition 3.1. Given a subexpression v of a reduced expression w = s il s i2 . . . s im , we define 



T° 


:={fc€{l,.. 


.,m} 


V( k - 


-i) 


< V(k)}, 


J? 


:={fcG{l,.. 


.,m} 




-i) 


= V(k)}, 




:={fc€{l,.. 


. ,ra} 


V( k - 


-i) 


> W(fc)}. 



The expression v is called non- decreasing if < V(j) for all j = 1, ... , to, e.g. J* = 0. 

The following definition is from [8, Definition 2.3]: 
Definition 3.2 (Distinguished subexpressions). A subexpression v of w is called distinguished if we have 

(3.1) v u) < s t] for all j e {1, . . . , m}. 

In other words, if right multiplication by Si j decreases the length of W(j-i), then in a distinguished 
subexpression we must have = v^-\)Si r 

We write v -< w if v is a distinguished subexpression of w. 

Definition 3.3 (Positive distinguished subexpressions). We call a subexpression v of w a positive dis- 
tinguished subexpression (or a PDS for short) if 

(3.2) v (j-i) < v (j-i) s ij for a11 J € {1, • • • ,m}. 
In other words, it is distinguished and non-decreasing. 

Lemma 3.4. [23] Given v < w and a reduced expression w for w, there is a unique PDS v + for v in w. 
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3.3. Deodhar components in the flag variety. We now describe the Deodhar decomposition of the 
flag variety. This is a further refinement of the decomposition of G/B into Richardson varieties 7Z V . W . 
Marsh and Rietsch [23] gave explicit parameterizations for each Deodhar component, identifying each 
one with a subset in the group. 

Definition 3.5. [23, Definition 5.1] Let w = s il 
distinguished subexpression. Define a subset G v 



(3.3) 



G v 



9i92 ■■■> 



There is an obvious map (M*)I J °I x 
v = w = 1 we define G v , w = {1}- 

Example 3.6. Let W = 65, w = S2S3S4S1S2S3 and v = s 2 llls 2 l. Then the corresponding element 
g £ G v , w is given by g = s 2 y 3 (p2)yi(p3)y 1 (pi)x 2 {m b )s^ 1 y d ,(p fi ) 1 which is 



x ie (me)s k 
VieiPe) 

Sit 



. Si m be a reduced expression for w, and let v be a 
in G by 
1 



if £e J*, 

if £e J°, 
if £e J°, 



for p< £ R*, ra f e R. 
,7 *l — >• G v , w defined by the parameters pi and me in (3.3). For 



/ 1 











o\ 


P3 


1 














Pa 


1 








P2P3 


p 2 - m 5 p e 


-m 5 


1 
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-P4P6 


-Pi 





V 



The following result from [23] gives an explicit paramctrization for the Deodhar component 7\L V 
We will take the description below as the definition of 7£ v ,w 



Proposition 3.7. [23, Proposition 5.2] The map 
isomorphism. The set G VjW lies in U~v<~)B + wB + , and the assignment g g-B + defines an isomorphism 



")I J °I x 



7v ' — > G v , w from Definition 3.5 is an 



(3.4) 

between the subset G v 



G v 



of the group, and the Deodhar component 1Z V w in G/B. 



Suppose that for each w £ W we choose a reduced expression w for w. Then it follows from Deodhar's 
work (see [8] and [23, Section 4.4]) that 



(3.5) 



7? 



□ ft v , w and G/B = \_\ I \_\ K v ^ j . 

v^w weW \v^w / 



These are called the Deodhar decompositions of1Z VtW and G/B. 

Remark 3.8. One may define the Richardson variety 1Z V . W over a finite field ¥ q . In this setting the 
number of points determine the i?-polynomials R v . w (q) = #(7?. 1 ,, t0 (F (? )) introduced by Kazhdan and 
Lusztig [14] to give a formula for the Kazhdan-Lusztig polynomials. This was the original motivation for 
Deodhar's work. Therefore the isomorphisms 7£ v ,w — (F* ) 1 x f|/ v ' together with the decomposition 
(3.5) give formulas for the i?-polynomials. 

Remark 3.9. Note that the Deodhar decomposition ofTZ VtW depends on the choice of reduced expression 
for w. However, we will show in Proposition 4.16 that its projection to the Grassmannian does not 
depend on the choice of reduced expression. 

Remark 3.10. The Deodhar decomposition of the complete flag variety is not a stratification - e.g. the 
closure of a component is not a union of components [10]. 

This decomposition has a beautiful restriction to the totally non-negative part (G/B)> of G/B. See 
[23, Section 11] and also [29] for more definitions and details. 
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Remark 3.11. Suppose we choose a reduced expression w for w, and for each v < w we let v + denote 
the unique positive distinguished subexpression for v in w. Note that v + is non-decreasing so J* = 0. 
Define G>° w to be the subset of G v+iW obtained by letting the parameters pi range over the positive 
reals. Let 7£>^ denote the image of G>° w under the isomorphism G v+ , w — ■ > 7^v+,w Then TJ^jJ, 
depends only on u and w, not on v + and w. Moreover, the totally non-negative part (G/B)> of G/B 
has a cell decomposition 



(3.6) (G/B)>q = U U ^ 

w£W \v<w 



>0 
v.w 



3.4. Projections of Deodhar components to the Grassmannian. The Richardson stratification of 
G/B has an analogue for partial flag varieties G/Pj introduced by Lusztig [22]. Let Wj be the parabolic 
subgroup of W corresponding to Pj, and let W J be the set of minimal- length coset representatives of 
W/Wj. Then for each w € W J , the projection tt : G/B — > G/Pj is an isomorphism on each Richardson 
variety 7\L„ iU ,. Lusztig showed that we have a decomposition of the partial flag variety 

(3.7) G/Pj = |_l ( U < n -^) 

w£W J \v<w 

Now consider the case that our partial flag variety is the Grassmannian Grk, n for k < n. The 
corresponding parabolic subgroup of W = & n is Wk = (si, «2, • • • , s n -k, ■ ■ ■ , s n -\). Let W k denote the 
set of minimal- length coset representatives of W/Wk- Recall that a descent of a permutation 7r is a 
position j such that > ir(j + 1). Then W k is the subset of permutations of & n which have at most 
one descent; and that descent must be in position n — k. 

Let TTk : G/B — > Grfe ; „ be the projection from the flag variety to the Grassmannian. For each 
w G W k and v < w, define V„ lW = ^k(J^v,w)- Then by (3.7) we have a decomposition 

(3.8) Gr Kn = U ( U V ^ 

wew k \v<w 

Remark 3.12. The decomposition in (3.8) coincides with the positroid stratification from Section 2.3. 
This was verified in [16, Theorem 5.9]. The appropriate bijection between the strata is defined in 
Lemma 3.13 below, and was first given in [34, Lemma A. 4]. 

Lemma 3.13. [34, Lemma A. 4] Let Q k denote the set of pairs (v, w) where v € W, w € W k , and v < w; 
let Dec„ denote the set of decorated permutations in S n with k weak excedances. We consider both 
sets as partially ordered sets, where the cover relation corresponds to containment of closures of the 
corresponding strata. Then there is an order-preserving bijection $ from Q k to Decjj which is defined as 
follows. Let (v,w) e Q J . Then &(v,w) — (n,col) where tt = vw^ 1 . We also let ir : (v,w) denote $>(v,w). 
To define col, we color any fixed point that occurs in one of the positions w(l), w(2), . . . , w(n — k) with 
the color —1, and color any other fixed point with the color 1. 

Since 7Tj. is an isomorphism from 1Z V . W to V VyW , it also makes sense to consider projections of Deodhar 
components in G/B to the Grassmannian. For each reduced decomposition w for w € W k , and each 
v -< w, we define T-V.w = TTfcC^-v.w)- Now if for each w G W k we choose a reduced decomposition w, 
then wc have 



(3.9) V v . w = |J "P v , w and Gr k , n 




10 



YUJI KODAMA AND LAUREN WILLIAMS 



Remark 3.14. By Remark 3.12 and Lemma 3.13, each projected Deodhar component P v 
positroid stratum S n -., where 7T = (it, col), it = uio -1 , and col is given by Lemma 3.13. 



lies in the 



Proposition 3.7 gives us a concrete way to think about the projected Deodhar components P v w . 
The projection ir k : G/B — > Gr kn maps each g e G v w to the span of its leftmost k columns. More 
specifically, it maps 

/Sl,n-fc+l ••• ffn,n-fe+l N 



' fi n 



Qn.n — k+1 



9 = 



9n,l 



fi.i, 



A = 



\ 9l,n 



9n.n 



\9l,n ■ ■ ■ 9l,n-k+l ■ 

Alternatively, we may identify A € Gr k ^ n with its image in the Pliicker embedding. Let a denote the 
column vector in R" such that the ith entry from the bottom contains a 1, and all other entries are 0, 
e.g. e„ = (1, 0, . . . , 0) T , the transpose of the row vector (1,0,..., 0). Then the projection ir k maps each 
g G G v , w (identified with g ■ B + e 7£ v , w ) to 



(3-10) 3 • e„ A . . . A e„_ fe+ i = £ Aj u ...,j k 

l<ji<...<jk<n 



{A)e jk A--- Ae h . 



That is, the Pliicker coordinate A, 



t (A) is given by 



■\7l, 



,3k * 



( A ) = ( e 3k A • • • A e ix , 3 • e„ A • • • A e„_ fe+ i) , 



where (•, •) is the usual inner product on 
Example 3.15. We continue Example 3.6. Note that w <G W k where k = 2. Then the map n 2 : G ViV 



Gr2,5 is given by 












f 1 











o\ 




P3 


1 











9 = 







1 










P2P3 




-m 5 


1 







V 




~P4 





V 



p 3 1 



4. Combinatorics of projected Deodhar components in the Grassmannian 

In this section we explain how to index the Deodhar components in the Grassmannian Grk. n by 
certain tableaux. We will display the tableaux in two equivalent ways - as fillings of Young diagrams 
by +'s and 0's, which we call Deodhar diagrams, and by fillings of Young diagrams by empty boxes, 
• 's and O's, which we call Go-diagrams. We refer to the symbols • and O as black and white stones. 

Recall that Wk = (s 1; s 2 , • • • , s n -k, ■ ■ ■ , s n-i) 1S a parabolic subgroup of W = & n and W h is the set 
of minimal- length coset representatives of W/Wk- 

An element w € W is fully commutative if every pair of reduced words for w are related by a sequence 
of relations of the form SiSj = SjSi. The following result is due to Stembridge [33] and Proctor [28]. 

Theorem 4.1. W k consists of fully commutative elements. Furthermore the Bruhat order on W k is a 
distributive lattice. 

Let Q k be the poset such that W k — J(Q k ), where J(P) denotes the distributive lattice of upper 
order ideals in P. The figure below (at the left) shows an example of the Young diagram of Gr^^. 
(The reader should temporarily ignore the labeling of boxes by Sj's.) The Young diagram should be 
interpreted as follows: each box represents an element of the poset Q k , and if b\ and 62 are two adjacent 
boxes such that 62 is immediately to the left or immediately above b\ , we have a cover relation b\ < 62 
in Q k . The partial order on Q k is the transitive closure of <. Note that the minimal and maximal 
elements of Q k are the lower right and upper left boxes, respectively. 
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We now state some facts about Q k which can be found in [33]. Let Wq € W k denote the longest 
element in W k . The simple generators used in a reduced expression for can be used to label Q k in 
a way which reflects the bijection between the minimal length coset representatives w £ W k and upper 
order ideals O w C Q k . Such a labeling is shown in the figure below. If b e O w is a box labelled by s i} 
we denote the simple generator labeling b by s;, := s,. Given this labeling, if O w is an upper order ideal 
in Q k , the set of linear extensions {e : O w — > [1, £(w)}} of O w are in bijection with the reduced words 
R(w) of w: the reduced word (written down from left to right) is obtained by reading the labels of O w 
in the order specified by e. We will call the linear extensions of O w reading orders. 



S5 


s 4 


S3 


S2 


Si 


S6 


S5 


S 4 


S3 


S2 


S7 


S6 


S5 


s 4 


S3 



15 


14 


13 


12 


11 


10 


9 


8 


7 


6 


5 


4 


3 


2 


1 



15 


12 


9 


6 


3 


14 


11 


8 


5 


2 


13 


10 


7 


4 


1 



Remark 4.2. The upper order ideals of Q k can be identified with the Young diagrams contained in 
a k x (n — k) rectangle, and the linear extensions of O w can be identified with the reverse standard 
tableaux of shape O w , i.e. entries decrease from left to right in rows and from top to bottom in columns. 

4.1. ©-diagrams and Deodhar diagrams. The goal of this section is to identify subexpressions of 
reduced words for elements of W k with certain fillings of the boxes of upper order ideals of Q k . In 
particular we will be concerned with distinguished subexpressions. 

Definition 4.3. [21, Definition 4.3] Let O w be an upper order ideal of Q k , where w £ W k . An ©-diagram 
("o-plus diagram") of shape O w is a filling of the boxes of O w with the symbols and +. 

Clearly there are 2*(w) ©-dia grams of shape O w . The value of an ©-diagram D at a box x is denoted 
D(x). Let e be a reading order for O w ; this gives rise to a reduced expression w = w e for w. The ffi- 
diagrams D of shape O w are in bijection with subexpressions v(D) of w: we will make the convention 
that if a box b <G O w is filled with a then the corresponding simple generator s b is present in the 
subexpression, while if b is filled with a + then we omit the corresponding simple generator. The 
subexpression v(D) in turn defines a Weyl group element v := v(D ) € W, where v < w. 

Example 4.4. Consider the upper order ideal O w which is Q k itself for 6 5 and k = 2. Then Q k is the 
poset shown in the left diagram. Let us choose the reading order (linear extension) indicated by the 
labeling shown in the right diagram. 



S3 


S2 


Si 


s 4 


S3 


S2 



6 


5 


4 


3 


2 


1 



Then the ©-diagrams given by 

























+ 











+ 






+ 





+ 





+ 



+ 


+ 





+ 





+ 



correspond to the expressions S2S3S4S1S2S3, IS3S4S1IS3, IS3IS1IS3, and IS3IS1II. The first and second 
are PDS's (so in particular are distinguished); the third one is not a PDS but it is distinguished; and 
the fourth is not distinguished. 
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Parts (1) and (2) of this proposition come from [21, Lemma 4.5 and Proposition 4.6]. 

Proposition 4.5. If b,b' G O w are two incomparable boxes, Sb and Sb> commute. Furthermore, if D is 
an 0-diagram, then 

(1) the element v :— v(D) is independent of the choice of reading word e. 

(2) whether v(D) is a PDS depends only on D (and not e). 

(3) whether v(D) is distinguished depends only on D (and not on e). 

Proof. The commutation of and st,> follows by inspection. For part (1), note that two linear 
extensions of the same poset (viewed as permutations of the elements of the poset) can be connected 
via transpositions of pairs of incomparable elements. Therefore v{D) is independent of the choice of 
reading word. 

Suppose D is an 0-diagram of shape O w , and consider the reduced expression w := w e = s il . . . Si n 
corresponding to a linear extension e. Suppose v(D) is a PDS of w. For part (2), it suffices to show 
that if we swap the fc-th and (k + l)-st letters of both w and v(Z)), where these positions correspond to 
incomparable boxes in O w , then the resulting subexpression v' will be a PDS of the resulting reduced 
expression w'. If we examine the four cases (based on whether the k-th and (k + l)-st letters of \r(D) 
are 1 or Si k ) it is clear from the definition that v' is a PDS. The same argument holds if v(D) is 
distinguished. □ 

This leads to the following definitions. Note that by Theorem 4.8, Definitions 4.6 and 2.13 agree. 

Definition 4.6. [21, Definition 4.7] A d-diagram of shape O w is an 0-diagram D of shape O w such that 
v(D) is a PDS. 

Definition 4.7. A Deodhar diagram of shape O w is an ©-diagram D of shape O w such that v(D) is 
distinguished. 

Theorem 4.8. [21, Theorem 5.1] and [27, Lemma 19.3] An 0-diagram is a J-diagram if and only if there 
is no which has a + above it (in the same row) and a + to its left (in the same column). 

Theorem 4.8 motivates the following open problem (which is slightly reformulated in Problem 4.13). 

Problem 4.9. Find an analogue of Theorem 4.8 for Deodhar diagrams which characterizes them by 
forbidden patterns. 

Definition 4.10. Let O w be an upper order ideal of Q k , where w € W and W — S n . Consider a 
Deodhar diagram D of shape O w ; this is contained in a k x (n — k) rectangle, and the shape O w gives 
rise to a lattice path from the northeast corner to the southwest corner of the rectangle. Label the steps 
of that lattice path from 1 to n; this gives a natural labeling to every row and column of the rectangle. 
We now let v be the permutation with reduced decomposition v(D), and we define n'-(D) to be the 
decorated permutation (tt(D), col) where n = n(D) = vw^ 1 . The fixed points of n correspond precisely 
to rows and columns of the rectangle with no +'s. If there are no +'s in the row (respectively, column) 
labeled by h, then n(h) = h and this fixed point gets colored with color 1 (respectively, — 1.) 

Remark 4.11. It follows from Remark 3.14 and the way we defined Deodhar diagrams that the projected 
Deodhar component Vd corresponding to D is contained in the positroid stratum S n -.m)- 

4.2. From Deodhar diagrams to Go-diagrams and labeled Go-diagrams. It will be useful for 
us to depict Deodhar diagrams in a slightly different way. Consider the distinguished subexpression v 
of w: for each k e J° we will place a O in the corresponding box; for each k € J* we will place a # in 
the corresponding box of O w ; and for each k € J° we will leave the corresponding box blank. We call 
the resulting diagram a Go-diagram, and refer to the symbols O and # as white and black stones. 
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Remark 4.12. Note that a Go-diagram has no black stones if and only if it corresponds to a Dcodhar 
diagram D such that v(D) is a PDS, i.e. a J-diagram. Therefore, slightly abusing terminology, we will 
often refer to a Go-diagram with no black stones as a J-diagram. 4 

Note that the Go-diagrams corresponding to the first three ©-diagrams in Example 4.4 are 



o 


o 


o 


o 


o 


o 



o 




o 


o 


o 





• 




o 




o 





Problem 4.13. Characterize the fillings of Young diagrams by blank boxes, white stones, and black 
stones which are Go-diagrams. 



Remark 4.14. Recall from Remark 3.8 that the isomorphisms 1Z V _ 



(F*)I J ° xf!,' 7 ' 1 together with the 
decomposition (3.5) give formulas for the i?-polynomials. Therefore a good combinatorial characteriza- 
tion of the Go-diagrams (equivalcntly, Deodhar diagrams) contained in a given Young diagram could 
lead to explicit formulas for the corresponding i?-polynomials. 

If we choose a reading order of O w , then we will also associate to a Go-diagram of shape O w a labeled 
Go-diagram, as defined below. Equivalcntly, a labeled Go-diagram is associated to a pair (v,w). 

Definition 4.15. Given a reading order of O w and a Go-diagram of shape O w , we obtain a labeled Go- 
diagram by replacing each O with a 1, each • with a —1, and putting a pi in each blank square b, where 
the subscript i corresponds to the label of b inherited from the linear extension. 

The labeled Go-diagrams corresponding to the examples above using the reading order from Example 
4.4 are: 



1 


1 


1 


1 — 1 


1 


1 



1 


P5 


1 


1 


1 


Pi 



-1 


P5 


1 


P3 


1 


Pi 



In future work we intend to explore further aspects of Go-diagrams and Deodhar strata. 

4.3. The projected Deodhar decomposition does not depend on the expressions w. Recall 
from Remark 3.9 that the Deodhar decomposition depends on the choices of reduced decompositions 
w of each w e W. However, its projection to the Grassmannian has a nicer behavior. 

Proposition 4.16. Let w G W k and choose a reduced expression w for w. Then the components of 
|J v ^ w 7^v,w do not depend on w, only on w. 

Proof. Recall from Theorem 4.1 that any two reduced expressions of w £ W k can be obtained from 
each other by a sequence of commuting moves (siSj — SjSi where \i — j\ > 2). And it is easy to check 
that if SiSj = SjSi, then 

(1) yi(a)yj(b) = yj{b)Vi{a) 

(2) 2/;(a)sj = Sjyi(a) 

(3) (xiia)^ 1 )^ = Sjixiia)^ 1 ) 

(4) (x i (a)sr 1 )y j (b)=y j (b)(x i (a)s^). 

The result now follows from Definition 3.5 and Proposition 3.7. □ 

4 Sincc J-diagrams arc a special case of Go-diagrams, one might also refer to them as Lego diagrams. 
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5. Plucker coordinates and positivity tests for projected Deodhar components 

Consider P ViW C Gr fei „, where w is a reduced expression for w G W k and v -< w. In this section 
we will provide some formulas for the Plucker coordinates of the elements of V VyVr , in terms of the 
parameters used to define G v , w - Some of these formulas are related to corresponding formulas for G/B 
in [23, Section 7]. 

5.1. Formulas for Plucker coordinates. 

Lemma 5.1. Choose any element A of ? v w c Gr kn . Then the lexicographically minimal and maximal 
nonzero Plucker coordinates of A are A/ and A// , where 

I = w {n, n — 1, . . . , n — k + 1} and I' = v {n, n — 1, . . . , n — k + 1}. 

Note that if we write / = . . . ,ik}, then /' = vw~ 1 {ii, . . . ,ik}- 

Proof. Recall that V VtV/ = 7r fc (7^ ViW ), where 7\L v .w C 1Z V _ W , and 1Z V . W = B + w ■ B + n B~v ■ B + . Now 
it is easy to check (and well-known) that the lexicographically minimal nonzero minor of each element 
in the Schubert cell iTk(B + w ■ B + ) is A/ and the lexicographically maximal minor of each element in 
the opposite Schubert cell irk(B~v ■ B + ) is A// where / and I' are as above. □ 

Our next goal is to provide formulas for the lexicographically minimal and maximal nonzero Plucker 
coordinates of the projected Deodhar components. 

Theorem 5.2. Let w — s il . . . Si m be a reduced expression for w e W k and v < w. Let / = w{n, n — 
1, . . . , n — k + 1} and I' — v{n, n — 1, . . . , n — k + 1}. Let A = nk(g) for any g e G v , w - If we write 
g = gi . . . g m as in Definition 3.5, then 

(5.1) Aj(A) = (-1)I J *I 11 Pl and A r (A) = 1. 

Note that A/ (A) equals the product of all the labels from the labeled Go-diagram associated to (v, w). 
Before proving Theorem 5.2, we record the following lemma, which can be easily verified. 

Lemma 5.3. For 1 < i < n — 1, we have 

(1) s^i = -e i+ i, s;e i+ i = e i; and s^j = ej if j ^ i or i + 1. 

(2) Ui(a)e l+ i = e i+ i + ae t and y l (a)e J = ej if j ^ i + 1. 

(3) (x l (a)s^ 1 )e i = e i+ \, (x i (a)s^ 1 )e i+ i = -(et + ae i+ i), and (x i (a)s~ 1 )e j = ej for j ^ i or i + 1. 
We now turn to the proof of Theorem 5.2. 

Proof. Recall from (3.10) how to identify each A <G Grk. n with its Plucker embedding. We first verify 
that A// (A) = 1. Since G v , w C U~i) (sec Proposition 3.7), we can write g G G v , w as g = hv with 
h e U~. Let A = e n A e„_i A • • • A e„_ fe+ i. Then A//(A) = (v ■ A, g • A) = (i) ■ A, hi) ■ A) = 1. 

Now we compute the value of A/(A). Recall from Proposition 4.16 that for w e W k , the Deodhar 
component 7\L v ,w does not depend on the choice of reduced expression w for w. Therefore we will fix a 
linear extension of Q k , and use that to construct our reduced expressions for each w € W k . 

It follows that each reduced expression w for w € W k where W = & n has the form 

(5.2) {Sj a S ja+1 . . . Sn-k+a-^iSja^Sj^+l . . . S n -fc+a-2) . . . {Sj 2 Sj 2+1 . . . S n _fe+i)(s jl S jl+1 . . . S„_fc). 

The four factors above correspond to the products of generators corresponding to the last, next-to-last, 
second, and top rows of the Young diagram, respectively. In particular, 1 < a < k (a is the number of 
rows in the Young diagram corresponding to w), and j\ < j% < ■ ■ ■ < j a -i < ja- Moreover, it is easy to 
check that {ji, J2, • ■ • , ja, n— k + a + l,n — k + a + 2, . . . ,n — 1, n} are the positions of the pivots of A 
(they correspond to the shape of the Young diagram) , so / = {ji, j 2 , ■ ■ ■ , j a ,n — k + a + l,n — k + a + 
2,. .. ,n - l,ra}. 
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Each g e G v , w will be obtained from (5.2) by replacing the Sj 



s by s,;'s, 



us write g = g^g^ ■ ■ .g^ where g^ is the product of g^s corresponding to {sj a Sj a+ i 
is the product of g^s corresponding to (sj a _ 



j/i(a)'s, or Xi(m)s i 's. Let 



-k+a-l) 



S ja-1+1 



n-fe+a-2)j etc. Now consider how such a g 
acts on e„, e n _i, .... Looking at Lemma 5.3, we see that g^ is the only portion of g which can affect 
e n -k+a (or any ej with j > n — k + a). This is because every Sj appearing in the other factors of (5.2) 
has the property that i < n — k + a — 2, and in this case, Sj, yi(a), and Xi(m)s~ l all act as the identity 
on e„_fc +a (or any ej with j > n — k + a). Similarly g^g^ is the only portion of g which can affect 
e n -fe+o-ij an d g^g^g^ is the only portion of g which can affect e„_fc +a ^ 2 , etc. 

Now we want to determine the value of the lexicographically minimal Plucker coordinate A/ (A). 
So we need to determine the coefficient of Ej in g ■ e n A ■ ■ ■ A e n -k+i- From Lemma 5.3, we see that 
s^ei+i = e,, ?/i(a)ei + i = ae^ + a higher term, and x i (a)s~' 1 ei + i = —e,i+ a higher term. Therefore from 
(5.2), we see that the expansion of g ■ e n -k+ a in the basis e\, . . . , e n has a nonzero coefficient in front of 
ej a . And that coefficient is (— l) q times the product of all the parameters p occurring in g^\ where q 
is the number of a;-factors in g^\ 

Similarly, from (5.2), the expansion of g ■ e n -k+a-i in the basis e\, . . . , e n has a nonzero coefficient 
in front of ej a _ 1 , and that coefficient is (— l) q times the product of all the parameters p occurring in 
g( 2 \ where q is the number of x-factors in g^ 2 \ 

Continuing in this fashion, the expansion of g ■ e n -k+i in the basis e\ , . . . , e„ has a nonzero coefficient 
in front of e^, and that coefficient is (— l) 9 times the product of all the parameters p occurring in g^ a \ 
where q is the number of x-factors in g^ a \ 

Additionally, g acts as the identity on e n -k+a+i, ■ ■ ■ , e„_i, and e„. It follows that the coefficient of 
Ej in the expansion of g • e n A • • • A e n -k+i in the standard basis is (— 1)' J «' Iliej Pii as desired. □ 

Our next goal is to give a formula for some other Plucker coordinates besides the lexicographically 
minimal and maximal ones. First it will be helpful to define some notation. 

Definition 5.4. Let W — & n , let w — s il . . . s im be a reduced expression for w e W k and choose v -< w. 
This determines a Go-diagram D in a Young diagram Y. Let b be any box of D. Note that the set 
of all boxes of D which are weakly southeast of b forms a Young diagram Y" b m ; also the complement of 
Y™ in Y is a Young diagram which we call F fe out (see Example 5.5 below). By looking at the restriction 
of w to the positions corresponding to boxes of Y" b m , we obtained a reduced expression w fc n for some 
permutation w b n , together with a distinguished subexpression v 6 n for some permutation v™. Similarly, 



by using the positions corresponding to boxes of Y b out , we obtained w£ ut , w b ut , v£ 



and v^ ut . 



When 



the box b is understood, we will often omit the subscript b. 

For any box b, note that it is always possible to choose a linear extension of O w which orders all the 
boxes of Y out after those of Y ln . We can then adjust w accordingly; Proposition 4.5 implies that this 
does not affect whether the corresponding expression v is distinguished. Having chosen such a linear 
extension, we can then write w = w m w out and v = v m v out . We then use g ln and <? out to denote the 
corresponding factors of g e G v , w . Wc define J^ out to be the subset of J° coming from the factors of 
v contained in v out . Similarly, for J° out and J* out . 

Example 5.5. Let W = 67 and w = S4S5S2S3S4S6S5S1S2S3S4 be a reduced expression for w e W 3 . 
Let v = S4S5IIS4IS5.S1IIS4 be a distinguished subexpression. So w = (3,5,6,7,1,2,4) and v = 
(2, 1, 3, 4, 6, 5, 7). We can represent this data by the poset O w and the corresponding Go-diagram: 



Si 


S3 


Si 


S\ 


Sb 


Si 


S3 


S2 


se 


S5 


Si 





• 






O 


• 


O 








O 


O 
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Let b be the box of the Young diagram which is in the second row and the second column (counting 
from left to right). Then the diagram below shows: the boxes of Y ln and Y out ; a linear extension which 
puts the boxes of Y out after those of Y' n ; and the corresponding labeled Go-diagram. Using this linear 
extension, w ln = S4S5S2S3S4, w out = SQS5S1S2S3S4, v m = S4S5IIS4, and v out = I.S5S1IIS4. 



out 


out 


out 


out 


out 


in 


in 


in 


out 


in 


in 





11 


10 


9 


8 


7 


5 


4 


3 


6 


2 


1 





-1 


P10 


P9 


1 


-1 


1 


P4 


P3 


P6 


1 


1 






Note that J* out = {7, 11} and J° out = {6, 9, 10}. Then g e G v , w has the form 

9 = 9 ln 9° nt = {sAS b y2{pz)yz{Pi)sA) (ye(P6)x 5 (m 7 )s b 1 s 1 y 2 (p9)y3(pio)x4(m 11 )s^ 1 ). 

When we project the resulting 7x7 matrix to its first three columns, we get the matrix 

-P3P10 -P10 -ma -1 0^ 
-P3P4 ~Pi ~m 7 1 
p 6 1, 

Theorem 5.6. Let w = . . . Si m be a reduced expression for w € W k and v -< w, and let D be the 
corresponding Go-diagram. Choose any box b of D, and let v ln = v b n and w m = w b n , and w out = v£ ut 
and w out = w% ut . Let A = n k (g) for any g e G v , w , and let / = w{n, n — 1, . . . , n — k + 1}. Define 
lb = v m (w m )^ 1 I G ('fe')- K we write g — gi . . . g m as in Definition 3.5, then 

(5.3) A 7b (A) = (-l)l^-tl [] 

Note that A/ b (A) equals the product of all the labels in the "out" boxes of the labeled Go-diagram. 

Remark 5.7. All of the Pliicker coordinates given by Theorem 5.6 are monomials in the p^s. In partic- 
ular, they are nonzero, and do not depend on the values of the m-parameters from the Xi (m)-factors. 

Before proving Theorem 5.6, we mention an immediate corollary. 

Corollary 5.8. Use the notation of Theorem 5.6. Let b be a box of the Go-diagram, and let e, s, and 
se denote the neighboring boxes which are at the east, south, and southeast of b. Then we have 

1 if box b contains a white stone 

— 1 if box b contains a black stone 

Pb if box b is blank and the labeled Go diagram contains Pb 

Remark 5.9. Each black and white stone corresponds to a two-term Pliicker relation, that is, a three- 
term Pliicker relation in which one term vanishes. And each black stone implies that there are two 
Pliicker coordinates with opposite signs. This will be useful when we discuss the regularity of solitons 
in Section 12. 

Example 5.10. We continue Example 5.5. By Theorem 5.2, / = w{5,6, 7} = {1,2,4} and /' = 
w{5,6, 7} = {5,6,7}, and the lexicographically minimal and maximal nonzero Pliicker coordinates for 
A are A/(A) = PsPaPqPsPio and Ai>(A) = 1; this can be verified for the matrix A above. 

We now verify Theorem 5.6 for the box b chosen earlier. Then lb = v ln (w ln )^ 1 I = {1,4, 5}. Theorem 
5.6 says that Ai b (A) = peP9Pio- The analogous computations for the boxes labeled 7, 6, 4, 3, 2, 
1, respectively, yield Ai i5i7 = -p 9 p w , Ai j2 ,7 = P3P4P9P10, Ai ;4i5 = p 6 p 9 pi , Ai i3;4 = p 4 P6P9Pio, 
Ai,2,4 = P3PaP&P9Pwi and A 1;2j 4 = P3PaP&P9Pw- These can be checked for the matrix A above. 



A Je (A)A Js (A) 
A Ib (A)&UA) 
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5.2. The proof of Theorem 5.6. For simplicity of notation, we assume that when we write A in 
row-echelon form, its first pivot is i\ = 1 and its last non-pivot is n. (The same proof works without 
this assumption, but the notation required would be more cumbersome.) 

Choose the box b which is the located at the northwest corner of the Young diagram obtained by 
removing the topmost row and the leftmost column; this is the box labeled 5 in the diagram from 
Example 5.5. We will explain the proof of the theorem for this box b. The same argument works if b 
lies in the top row or leftmost column; and such an argument can be iterated to prove Theorem 5.6 for 
boxes which are (weakly) southeast of b. 

Choose a linear extension of O w which orders all the boxes of Y out after those of Y in , and which orders 
the boxes of the top row so that they come after those of the leftmost column. The linear extension 
from Example 5.5 is one such an example. Choosing the reduced expression w correspondingly, we 
write w = w ln w out and v = v ln v out , then choose g G G v , w an d write it as g = g ln g out . Note that from 
our choice of linear extension, we have 

(5.4) w out = (s„_is„_ 2 • • • s„_ fe+ i)(sis 2 ■ • ■ s n -k)- 

Recall that Ib — v ln (w ln )~ 1 I where / = . . . , i^}, with i\ = 1. Also iu _1 J = {n— fc+1, . . . , n— 1, n} 7 
which implies that 

(5.5) (tu™)- 1 / = w out {n - k + 1, . . . , n - 1, n} = {1, n - k + 1, n - k + 2, . . . , n - 1}. 

Since there is no factor of si or s„_i in v m , and = w m {l, n — k + 1, n — k + 2, . . . , n — 1}, we have 

(5.6) 1 e h and n ^ I b . 

Write I b = {ji, . . . , j fe } with ji = 1. Our goal is to compute A Ib (A) = (e jk A- ■ ■f\e jl ,g-e n f\- ■ -Ae„_ fe+ i). 

Let fi — g ■ e n -k+i- Let qi be the product of all labels in the "out" boxes of the £th row of the 
labeled Go-diagram. Using Lemma 5.3 and equation (5.4), we obtain 

fk = 9 ■ e„ = g in ■ (qke n -i + c k n e n ) 
fk-i = 9 ■ e n -i = 9 m ■ (9fc-ie n _ 2 + c^lie„_i + c£ -1 e„) 

h = 9 ■ en-k+2 = 9 m ■ (<72e„-fc+i + c 2 n _ k+2 e n -k+2 H h c 2 n e n ) 

h = 9 ■ e n -k+i = 9 m ■ (<7iei + c\e 2 H h c\e n ). 

We now claim that only the first term with coefficient qi in each fe contributes to the Pliicker 
coordinate Af,(^4). To prove this claim, note that: 

(1) Since n ^ I h and g ln ■ e n = e„, the terms c l n e n do not affect A b (A). Therefore, we may as well 
assume that each c l n — 0. Define fk — qk9 m 1 e„_i. 

(2) Now note that the term c n Z\e n -i docs not affect the wedge product fk A fk-i- In particular, 
fk A fk-i = fk A fk-i where f k -i = qk-ig ln ■ e n -2- 

(3) Applying the same argument for 2 < I < k — 2, we can replace each ft by f e = q e g m ■ e n -k+£, 
without affecting the wedge product. 

(4) Since 1 <E h and e\ does not appear in any fi except /i, for the purpose of computing Ab(^4) 
we may replace f Y by f 1 = q x e x - 
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Now we have 

&h( A ) = ( e « A-Ae^/nA-A/i) 
= (ej k A---Ae jl J k A---Af 1 ) 

k 

= (11*)^ A • • • A e jn9 ln ■ (e„-i A • • • A e„_ fe+1 A d)) 
j=i 
fe 

= (11*)^ A "' Ae J2,9 ln ■ (e n -i A •• • Ae„_ fe+ i)), 

where in the last step we used ji = 1. But now note that 

(e jk A • • • A e h ,g m • (e„_i A • • • A e„_ fe+1 )) = 1, 

because this is the lexicographically maximal minor for the matrix Ab = irk-i(g ln ) € Gr k -i , n -2 corre- 
sponding to the sub Go-diagram obtained by removing the top row and leftmost column. 
Therefore A h (A) = n} =1 Qj = (-l) |7 * out 1 Etiej^ Pi, as desired. 

5.3. Positivity tests for projected Deodhar components in the Grassmannian. We can use our 

results on Pliicker coordinates to obtain positivity tests for Deodhar components in the Grassmannian. 

Definition 5.11. Let D be a Go-diagram and Sd C Gr k ^ n . A collection J of fc-clcmcnt subsets of 
{1,2,..., n} is called a positivity test for So if for any A G Sx>, the condition that Ai(A) > for all 
I & J implies that A e (GVfc ; „)>o- 

Theorem 5.12. Consider A G Grk tU lying in some Deodhar component 5d, where D is a Go-diagram. 
Consider the collection of minors J = {A](A)}U{A] b (A) | b a box of D}, where I and lb are defined as 
in Theorem 5.6. If all of these minors are positive, then D has no black stones, and all of the parameters 
Pi must be positive. It follows that the Deodhar diagram corresponding to D is a J-diagram, and A 
lies in the positroid cell 5^" c {Gr k ,„)>o- In particular, J' is a positivity test for Sp. 

Proof. By Remark 5.9, if all the minors in J are positive, then D cannot have a black stone. 
By Theorem 5.2 and Theorem 5.6 we have that 

Aj(A) = (-1)I J *I Y[ Pi and A h (A) = (-1)1^1 



„out 



Since we are assuming that both of these are positive, it follows that for any box b, we have that 

is also positive. Now by considering the boxes b of D in an order proceeding from southeast to northwest, 
it is clear that every parameter pi in the labeled Go-diagram must be positive, because each ^ 
must be positive. 

Let v and w be the Weyl group elements corresponding to D. Then it follows from Remark 3.11 
that A lies in the projection of the totally positive cell 1Z>° W . And the projection of TZ^° W is precisely 
the positroid cell S^™ of (Grk, n )>o- □ 

6. SOLITON SOLUTIONS TO THE KP EQUATION 

We now explain how to obtain a soliton solution to the KP equation from a point of Gr k , n - Each 
soliton solution can be considered as an orbit with the flow parameters (x,y,t) € R 3 on Gr fci „. 



THE REGULARITY OF KP SOLITONS 



19 



6.1. Prom a point of the Grassmannian to a r-function. We start by fixing real parameters Kj 
such that 



K\ < n 2 < 



< K r , 



which are generic, in the sense that the sums K j m are an distinct for any p with 1 < p < n. We 

m— 1 

also assume that the differences between consecutive k,'s are similar, that is, — is of order one 
(e.g. one can take all Kj to be integers). 

We now give a realization of GV^„ with a specific basis of E". We define a set of vectors {E° : j = 
1,. .. ,n} by 



E u := 



.n-2 



V 1 / 

Since all Kj-'s are distinct, the set {E° : j = 1, . . . , n} forms a basis of R n . Now define an n x n matrix 
E° = (E° , . . . , E°), and let A be a full-rank k x n matrix parametrizing a point on Grk, n - Then the 
vectors {F° e M™ : i = 1, . . . , k} span a /c-dimensional subspace in R", where F° is defined by 



F° — V a ■ E° 

3 = 1 



or (F?,...,F°)^i? ^. 



For I = {ii, . . . , i k ), define the vector E° = E° A • • • A E° . Then we have a realization of the Pliicker 
embedding: 

F?A---AF°= E Aj(A)E?. 

In [30] , Sato showed that each solution of the KP equation is given by an orbit on the Grassmannian. 
To construct such an orbit, we consider a deformation E* of the vector E°, defined by: 

t:=(x,y,z), 8j(x,y,z) = Kjx + rfy + Kjt, E* := E° exp (6j(x, y, tj) . 

Remark 6.1. Let E t be the n x n matrix function whose columns are the vectors {E*}: 

£7* := (E* l7 . . . , E*) = J B°diag(e e ", e e "-, . . . ,e e >). 

Note that E° is a Vandcrmonde matrix. The vector functions {E*} form a fundamental set of solutions 
of a system of differential equations. More concretely, if we define elementary symmetric polynomials 
in the Kj 's by 



Ol = J2 K 3> a 2 — J2 K i K 3> 
J=l i<j 



&3 



y~] KiKjKk, 
i<j<k 



and let Ck be the companion matrix 



a 



K 





-cr 2 • • 


• • • • ±CT„\ 


l 





• • •• 





1 


: 






• : 







• 1 J 



then the matrix E t satisfies 



CE t := 



dx 



— Ck \ E = 0. 
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So for any t = (a;, y, t), we have 

E" £* kcr(£) = Span R {E* : j = 1, . . . , n}. 

Note that can be diagonalized by the Vandermonde matrix E°, i.e. 

C K £° = £°A where D = diag(«„, . . . , m). 

Each vector function E t satisfies the following linear equations with respect to y and t: 
dE t d 2 E t QE t Q3 E t 

-W = ^ = CkE and ~W ~ ~ CkE ■ 

This is a key of the "integrability" of the KP equation, that is, the solutions of the linear equations 
provide a solution of the KP equation. 

We now define an orbit generated by the matrix E t on elements of G — GL„, 

g t := E t g for each g g GL„ . 

Then {g* • e„ A • • • A e n -k+i} is a flow (orbit) of the highest weight vector on the fundamental repre- 
sentation space of GL„. 

Next we define the r-function as 

r(x,y,t) := (e k A •• • Aei, g t ■ e n A •• • A e n _ fe+ i). 

Given I = {ii, . . . , ik} C [n], we let £J/(x, y, i) denote the scalar function 

E^x^t) =<e fc A--- Aei, A-AEjJ =<e fc A---Aei, E? A • • • A E9 ) e o*i +-+««* 

(fU) =fll(^-^))^ + - +e - fc - 

With the projection 7Tfe : SL„ — > Grk. n , g ^ A, the r-function can be also written as 
(6.2) r(a:,i/,t)=TA(a;,tf,t)= E Aj(^)£?j(a;,y,t). 

It follows that if A g (Gr fe ,„)> , then t a > for all (x,y,t) g M 3 . 

Remark 6.2. The present definition of the r-function is quite useful for the study of the Toda lattice 
whose solutions are defined on a complete flag manifold. We will discuss the totally non-negative flag 
variety and the Toda lattice in a forthcoming paper. 

6.2. Prom the r-function to solutions of the KP equation. The KP equation for u(x,y,t) 

d ( ^du du d 3 u\ 2,® 2u 
dx \ dt dx dx 3 J dy 2 

was proposed by Kadomtsev and Petviashvili in 1970 [13], in order to study the stability of the soliton 
solutions of the Korteweg-de Vries (KdV) equation under the influence of weak transverse perturbations. 
The KP equation can be also used to describe two-dimensional shallow water wave phenomena (see for 
example [18]). This equation is now considered to be a prototype of an integrable nonlinear partial 
differential equation. For more background, see [25, 9, 1, 12, 24]. 

Note that the r-function defined in (6.2) can be also written in the Wronskian form 

(6-3) T A (x,y,t)=W(f u f 2 ,...,f k ), 
with the scalar functions {fj : j = 1, . . . , k} given by 

(/i,/2,- • • ,.fk) T = A ■ (exp0i,exp0 2 , ■ • • ,exp6»„) T , 
where (. . .) T denotes the transpose of the (row) vector (...). 
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It is then well known (see [12, 5, 6, 7]) that for each choice of constants {ki, . . . , K n } and element 
A £ Grk,m the r-function defined in (6.3) provides a soliton solution of the KP equation, 

d 2 

(6.4) u A (x,y,t) = 2 — In T A (x,y,t). 

If A £ (Grk,n)>0) then it is obvious that u A (x,y,t) is regular for all (x,y,t) £ K 3 . A main result of 
this paper is that the converse also holds - see Theorem 12.1. Throughout this paper when we speak 
of a soliton solution to the KP equation, we will mean a solution u A (x,y, t) which has the form (6.4), 
where the r-function is given by (6.2). 

Remark 6.3. The function Ei(x,y,t) in the r-function (6.2) can be expressed as the Wronskian form 
in terms of {Ei j = e 9 ^ : j = 1, . . . , fc}, i.e. 

E!(x,y,t) =Wv(E h ,E i2 ,. . . ,E ih ). 



7. Contour plots of soliton solutions 

One can visualize a solution u A (x,y,t) to the KP equation by drawing level sets of the solution in 
the rry-plane, when the coordinate t is fixed. For each r £ R, we denote the corresponding level set by 

C r (t) := {(x,y) £ M 2 : u A (x,y,t) = r}. 

Figure 2 depicts both a three-dimensional image of a solution u A (x, y,t), as well as multiple level sets 
C r . These level sets are lines parallel to the line of the wave peak. 





FIGURE 2. A line-soliton solution u A (x,y,t) where A — (1,1) £ (Gri^ 2 )>o, depicted 
via the 3-dimensional profile u A (x, y, t), and the level sets of u A (x,y,t) for some t. Ei 
represents the dominant exponential in each region. 



To study the behavior of u A (x,y,t) for A £ Sm C Gr^.m we consider the dominant exponentials in 
the r-function (6.2) at each point (x,y,t). First we write the r-function in the form 

r A (x,y,t)= £ Aj{A)Ej(x,y,t) 



= E exp E {K h x + K 2 n y + Kp) + MAj(A)Kj) 

J£M 

where Kj := Y\i <m { K j m ~ K j e ) > 0. Note that in general the terms ln(Aj(A)Kj) could be imaginary 
when some Aj(A) are negative. 
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Since we are interested in the behavior of the soliton solutions when the variables (x,y,t) are on a 
large scale, we rescale the variables with a small positive number e, 



x y t 

y — > -, t — ► -. 
e e e 



This leads to 



rl(a;,2/,*)= £ exp ( - £ ( Kji x + K%y + K%t) + \n(Aj(A)K, ; ) ) . 

JEM \ e »=1 J 

Then we define a function Ja(x, y, t) as the limit 

f A {x,y,t) = ]imeln (T%(x,y,t)) 

f£, 2 3 ^ 

= max | E («i 4 a; + 4 2/ + • 

Since the above function depends only on the collection A4, we also denote it as fM{x,y,t). 

Definition 7.1. Given a solution UA{x,y,t) of the KP equation as in (6.4), we define its contour plot 
C(ua) to be the locus in R 3 where fA(x,y,t) is not linear. If we fixt = to, then we let Ct (uA) be the 
locus in M 2 where f A {x, y,t = to) is not linear, and we also refer to this as a contour plot. Because these 
contour plots depend only on M. and not on A, we also refer to them as C(Ai) and C to (Ai). 

Remark 7.2. The contour plot approximates the locus where \uA(x,y,t)\ takes on its maximum values 
or is singular. 

Remark 7.3. Note that the contour plot generated by the function f A {x, y,t) at t = consists of a set 
of semi-infinite lines attached to the origin (0, 0) in the a;j/-plane. And if t\ and ti have the same sign, 
then the corresponding contour plots C tl (A4) and C t2 (M) are self-similar. 

Also note that because our definition of the contour plot ignores the constant terms ln(Aj(A)Kj), 
there are no phase-shifts in our picture, and the contour plot for fA(x,y,t) = fM( x ,y,t) does not 
depend on the signs of the Pliicker coordinates. 

It follows from Definition 7.1 that C(ua) and C to (uA) are piecewise linear subsets of M 3 and M 2 , 
respectively, of codimension 1. In fact it is easy to verify the following. 

Proposition 7.4. [20, Proposition 4.3] If each Ki is an integer, then C{ua) is a tropical hypersurface in 
M 3 , and Ct Q {uA) is a tropical hypersurface (i.e. a tropical curve) in R 2 . 

The contour plot Ct (uA) consists of line segments called line-solitons, some of which have finite 
length, while others are unbounded and extend in the y direction to ±oo. Each region of the complement 
of Ct (ua) in M 2 is a domain of linearity for Ja{x, y, t), and hence each region is naturally associated to 
a dominant exponential Aj(A)Ej(x,y,t) from the r-function (6.2). We label this region by J or Ej. 
Each line-soliton represents a balance between two dominant exponentials in the r-function. 

Because of the genericity of the K-parameters, the following lemma is immediate. 

Lemma 7.5. [7, Proposition 5] The index sets of the dominant exponentials of the r-function in adjacent 
regions of the contour plot in the xy-plane are of the form {i, l 2 , ■ ■ ■ , Ik} and {j, I2, ■ ■ ■ , h}- 

We call the line-soliton separating the two dominant exponentials in Lemma 7.5 a line-soliton of type 
Its equation is 

(7.2) x + (Ki + Hj)y + (k\ + + K 2 )t = 0. 

Remark 7.6. Consider a line-soliton given by (7.2). Compute the angle between the positive y-axis 
and the line-soliton of type measured in the counterclockwise direction, so that the negative x-axis 
has an angle of ^ and the positive x-axis has an angle of — |. Then tan^^ jj = Ki + kj, so we refer to 
^ + Kj as the slope of the line-soliton (see Figure 2). 
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In Section 9 we will explore the combinatorial structure of contour plots, that is, the ways in which 
line-solitons may interact. Generically we expect a point at which several line-solitons meet to have 
degree 3; we regard such a point as a trivalent vertex. Three line-solitons meeting at a trivalent vertex 
exhibit a resonant interaction (this corresponds to the balancing condition for a tropical curve). See 
[20, Section 4.2]. One may also have two line-solitons which cross over each other, forming an X-shape: 
we call this an X-crossing, but do not regard it as a vertex. See Figure 4. Vertices of degree greater 
than 4 are also possible. 

Definition 7.7. Let i < j < k < I be positive integers. An X-crossing involving two line-solitons of 
types [i, k] and [j, £] is called a black X-crossing. An A-crossing involving two line-solitons of types [i, j] 
and [k,£], or of types [i,£] and [j, fc], is called a white X-crossing. 

Definition 7.8. A contour plot C t {uA) is called generic if all interactions of line-solitons are at trivalent 
vertices or are A-crossings. 

Example 7.9. Consider some A e Gr 4i g which is the projection of an element g <G G v , w with 

w = s 7 ssS4S5SQSrS2S4,S5SeSiS2S 3 S4S 5 and v = S7IIS5IS7S2IS4IHS2IS4S5. 
Then v = 1 and n = vw^ 1 = (6, 7, 1, 8, 2, 3, 9, 4, 5). The matrix g e G v , w is given by 

.9 = S7j/8(P2)j/4(P3)S52/6(P5)aJ7(m 6 )sf 1 S22/3(P8)S4?/5(Pl0)j/6(Pll) 

• yi(pi2)x 2 (m 13 )s2 1 y 3 (p 14 )x 4 (m 15 )s4 1 x 5 (m 16 )s^ 1 . 
The Go-diagram and the labeled Go-diagram are as follows: 



• 


• 




• 
















• 























-1 


-1 


Pl4 


-1 


Pl2 


Pll 


P10 


1 


P8 


1 


-1 


P5 


1 


P3 




P2 


1 









The A-matrix is then given by 
/ 

A = 



-P12P14 






913 

P8P10P11 





Pli 






915 

Pll(.P3 + Pw) 
-P3P5 





-77116 
Pll 
~P5 









where the matrix entry (713 = —mispu 





V 

f pro). In Figure 3, 



1 
1 

— 777g 
P2 

- wisps ~ m le p$pi and 1715 = m 15 - m le (p 3 
we show contour plots Cj(ua) for the solution 77,4(2^ y, t) at t = —10, 0, 10, using the choice of parameters 
. . . , k 9 ) = (—5, —3, —2, —1, 0, 1, 2, 3, 4), pj = 1 for all j, and to/ = for all t. Note that: 

(a) For y ^> 0, there are four unbounded line-solitons, whose types from right to left are: 

[1,6], [2,7], [4,8], and [7,9] 

(b) For y <C 0, there are five unbounded line-solitons, whose types from left to right are: 

[1,3], [2,5], [3,6], [4,8], and [5,9] 

Apparently the line-solitons for y >• correspond to the excedances in tt = (6, 7, 1, 8, 2, 3, 9, 4, 5), while 
those for y <C correspond to the nonexcedances. In Section 8 we will give a theorem explaining the 
relationship between the unbounded line-solitons of Ct(uA) and the positroid stratum containing A. 

Note that if there are two adjacent regions of the contour plot whose Plucker coordinates have different 
signs, then the line-soliton separating them is singular. For example, the line-soliton of type [4,8] (the 
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t = - 10 t = t= 10 



Figure 3. Example of contour plots Ct{uA) for A £ Gr^g. The contour plots are 
obtained by "Plot3D" of Mathematica (see the details in the text). 

second soliton from the left in y 3> 0) is singular, because the Pliicker coordinates corresponding to the 
(dominant exponentials of the) adjacent regions are 

Ai, 2 ,4,9 = P3P5PsPiaPnPi2Pi4 = 1 and Ai i2)8 ,9 = -P8P10P11P12PU = -1. 

8. Unbounded line-solitons at y > and y <c 

In this section we show that the unbounded line-solitons at \y\ 3> of a contour plot Ct{uA) are 
determined by which positroid stratum contains A. Conversely, the unbounded line-solitons of Ct{uA) 
determine which positroid stratum A lies in. The main result of this section is Theorem 8.1. 

Theorem 8.1. Let A £ Gr kn lie in the positroid stratum S„: , where tt — (7r, col). Consider the contour 
plot Ct{uA) for any time t. Then the excedances (respectively, nonexcedances) of 7r are in bijection with 
the unbounded line-solitons of C t {uA) at y ^> (respectively, y <C 0). More specifically, in C t (uA), 

(a) there is an unbounded line-soliton of [i, /i]-type at y 3> if and only if n(i) = h for i < h, 

(b) there is an unbounded line-soliton of [i, /i]-type at y <C if and only if 7r(/i) = i for i < h. 

Therefore w determines the unbounded line-solitons at y ^> and y <C of Ct(uA) for any time t. 

Conversely, given a contour plot Ct(uA) at any time t where A € Gr^n, one can construct w = (it, col) 
such that A £ S^-. as follows. The excedances and nonexcedances of tt are constructed as above from 
the unbounded line-solitons. If there is an h £ [n] such that h £ J for every dominant exponential Ej 
labeling the contour plot, then set tt(Ii) = h with col(h) = 1. If there is an h £ [n] such that h ^ J for 
any dominant exponential Ej labeling the contour plot, then set ir(h) = h with col(h) = —1. 

Proof. This result will follow immediately from Theorems 8.3 and 8.7 below. □ 

Remark 8.2. Chakravarty and Kodama [5, Prop. 2.6 and 2.9] and [7, Theorem 5] associated a derange- 
ment to each irreducible element A in the totally non-negative part (Grk, n )>o of the Grassmannian. 
Theorem 8.1 generalizes their result by dropping the hypothesis of irreducibility and extending the 
setting from (Gr k , n )>Q to Gr k , n . 

Before stating Theorems 8.3 and 8.7, we need to introduce some notation. 

Given a matrix A with n columns, let A(k, ...,£) be the submatrix of A obtained from columns 
k, + . . . , I— l,t, where the columns are listed in the circular order k, fc + 1, . . . ,n— 1, n, 1, 2, . . . , k — 1. 

The following result generalizes [2, Lemma 3.4] from (Gri~,n)>o to Gr k n . Our proof of Theorem 8.3 
will be similar to that of [2], but some arguments can be clarified using some basic theory of matroids. 
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Theorem 8.3. Let A G Gr^^ n and consider the contour plot Ct{uA) for any time t. Choose i,h € 
{1, . . . , n} with i < h. 

Then there is an unbounded line-soliton of Ct(v,A) at y labeled [i, h] if and only if 

(8.1) rank A(i, . . . ,h - 1) = rank A(i + 1, . . . , h) = rank A(i, ...,h)= rank A(i + 1, . . . , h - 1) + 1. 
And there is an unbounded line-soliton of Ct{uA) at y ^> labeled [i, ft-] if and only if 

(8.2) rank A(h, ...,i-l) = rank A(h + l,...,i) = rank A(h, ...,i) = rank A(h + 1, . . . , i - 1) + 1. 

Recall from Section 6 that 0j(ar, y, z) = KjX + K^y + Kjt. Fix i, j S {1, . . . , n}, and let Ljj denote the 
line defined by 0j = Oj. Define subsets of [n] by 

P = {max(i, j) + 1, . . . , niin(i,j) — 1} := {1,..., mm(i,j) — 1} U {max(«, j) + 1, . . . ,n} and 

Q = {min(i, j) + 1, . . . ,max(i, j) - 1}. 

In order to study the unbounded solitons at y 3> and y <§C 0, we first record the following lemma. 

Lemma 8.4. [2, Lemma 3.1] For |y| 3> 0, we have the following ordering among the 0j's on the line Lij: 

(1) For y 3> on the line Lij, 9 m < 8i = 9j for all m e Q, and 6* m > 0, = 0j for all m e P. 

(2) For y <C on the line Ljj, m > ^ = 6*^ for all m e Q, and 6* m <@i = 9j for all m e P. 

Proof. For a fixed t, the equation of the line Ljj (which is defined by ^ = Oj) has the form 

x + («i + Kj)y = constant. 

Then along Lj^, we have 

9m — ^m' = ( K m — K m')[( K m "F K rn') — ( K i + K j)]?/ + S, 

where 6 does not depend on x or y. The lemma now follows from the fact that K\ < K2 < ■ ■ ■ < K n . □ 
Then it follows immediately that 

Corollary 8.5. For y = y 3> (respectively y = yo <C 0) there is a well-defined total order on 6\, . . . , 6 n 
on the line Lij (with 9i — 9j), and this order does not change if we increase y (resp., decrease y). 

The following matroidal result will be useful to us. 

Proposition 8.6. [26, Theorem 1.8.5] Consider a matroid M. of rank k on the set [n], and let oj = 
(oji, . . . ,u> n ) € R n . Define the weight of a basis J = (ji, • • ■ ,jk) of M. to be u>j 1 + • • • + u>j k . Then the 
basis (or bases) of maximal weight are precisely the possible outcomes of the greedy algorithm: Start 
with J = 0. At each stage, look for an w-maximum element of [n] which can be added to J without 
making it dependent, and add it. After k steps, output the basis J. 

We now turn to the proof of Theorem 8.3. We will prove the result for unbounded line-solitons at 
i/>0 (the other part of the proof is analogous). 

Proof. Let M. be the matroid associated to A. Its ground set [n] is identified with the columns of A. 
First suppose that for i,j G [n], with i > j we have 

(8.3) rank A(i, . . . , j - 1) = rank^(i + l,...,j) = rank A(i, ...,j) = rank A(i + 1, . . . , j - 1) + 1. 

By Corollary 8.5, at y > we have a well-defined total order on the 8 m } s on the line L iy j. At y 3> 
the problem of computing the dominant exponential is equivalent to finding the basis of M with the 
maximal weight with respect to (9\, . . . , 6 n ). 

By Proposition 8.6, we can compute such a weight-maximal basis using the greedy algorithm. By 
Corollary 8.4, the greedy algorithm will first choose as many columns of A(i+1, . . . , j — 1) as possible. All 
of the # m 's are distinct except for 9i = 6j, so there will be a unique way to add a maximal independent 
set of columns of A(i + 1, . . . , j — 1) to the basis we are building. Note that by (8.3), the rank of 
A(i + 1, . . . , j — 1) is less than k, so our weight-maximal basis must additionally contain at least one 
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column that is not from A(i + 1, . . . , j — 1). By Corollary 8.4, columns i and j share a weight which is 
greater than any of the other remaining columns, so the next step is to add one of columns i and j to 
the basis we are building. By (8.3), we cannot add both columns, because doing so will only increase 
the rank by 1. Therefore we now have two ways to build a weight-maximal basis, by adding either one 
of the columns i and j. If the two bases we are building do not yet have rank k, then there is now a 
unique way to add columns from A(j + 1, . . . , i — 1) to complete both of them. 

We have now shown that along L^j at y 3> 0, there are precisely two dominant exponentials, Ej and 
Ej, where / = ( J U {i}) \ {j}. Therefore there is an unbounded line-soliton at y ^> labeled [j, i\. 

Conversely, suppose that for i > j, there is an unbounded line-soliton labeled [j, i] at y ^> 0. Then on 
the line L iy j there are two dominant exponentials Ej and Ej with J = (I Li {j}) \ {i}. By Proposition 
8.6, these must be the two outcomes of the greedy algorithm. As before, by Corollary 8.4, the greedy 
algorithm will first choose as many columns of A(i + 1, . . . , j — 1) as possible while keeping the collection 
linearly independent, and then the next step will be to add exactly one of the columns i and j. Since 
neither dominant exponential contains both i and j, adding both columns must not increase the rank 
more than adding just one of them. Therefore equation (8.3) must hold. □ 

Theorem 8.7. Let A € Gr^^n he in the positroid stratum S n -. where 7T = (it, col). Choose 1 < i < h < n. 
Then Tr(h) = i if and only if equation (8.1) holds, and n(i) — h if and only if equation (8.2) holds. 

Proof. Let I = (I\, . . . ,I n ) be the Grassmann necklace associated to A, so + = + {1). Then by 
Lemma 2.8, Ii — {xi,x 2 , . . . ,Xk} is the lexicographically minimal fc-subset with respect to the order 
i < i + 1 < ••• < n < 1 < ••• < i — 1 such that A/. (A) ^ 0. Similarly is the lexicographically 
minimal /c-subset with respect to the order i + 1 < ■ ■ ■ < n < 1 < ■ ■ ■ < i — 1 < i such that A/. +1 (^4) ^ 0. 

We will prove the first statement of the theorem (the proof of the second is analogous, so we omit 
it.) Suppose that Tr(h) = i. Then x\ — i; otherwise the ith column of A is the zero-vector and Tr(i) = i. 
Using Definition 2.14 and Lemma 2.8, h has the following characterization. Consider the column indices 
in the order i + 1, i + 2, . . . , n, 1, 2, . . . , i and greedily choose the earliest index h such that the columns 
of A indexed by the set {x 2 , . . . , Xk} U {h} are linearly independent. Then = (Ii \ {i}) U {h}. 

Now consider the ranks of various submatriccs of A obtained by selecting certain columns. 

Claim 0. rank A(i + 1, . . . , h — 1, h) = 1 + rank^4(i + 1, . . . , h — 1). This claim follows from the 
characterization of h and the fact that Jj + i is the lexicographically minimal /c-subset with respect to 
the order i + l<---<n<l<---<i such that A/. +1 (A) ^ 0. 

Claim 1. rank A(i, ... ,h) = rank A(i, h—1). To prove this claim, we consider two cases. 

Either i <i h <i Xk or i <i Xk <i h, where <i is the total order i<i + l<---<n<l<---<i — 1. 
In the first case, the claim follows, because h is not contained in the set Ii. In the second case, 
rank A(i, i+ 1, i + 2, . . . ,Xk) = k, and the index set {i, i + 1, . . . ,Xk} is a proper subset of {i, i + 1, . . . , h}, 
so rank A(i, . . . , h) = rank A(i, . . . , h — 1) = k. 

Now let R = rankA(i + 1, i + 2, . . . , h - 1). By Claim 0, rank A(i + 1, . . . , h) = R + 1. Therefore we 
have rank A(i, . . . , h) > rankA(z + l, . . . , h) = R+l. By Claim 1, rank A(i, . . . , h) — rank^4(i, . . . , h—1), 
but rankyl(i, . . . , h — 1) < R + 1, so rank^4(i, . . . , h) < R + 1. We now have rankA(i, . . . , h) = R+l. 
But also rank A(i, . . . , h — 1) = rank A(i, . . . , h) = R+l. Therefore rank^4(i, i + 1, . . . , h — 1) = 
rank A(i + 1, . . . , h — 1, h) = rank A(i, . . . , h) = rank^4(i + 1, . . . , h — 1) + 1, as desired. 

Conversely, suppose that rank A(i, i + 1, . . . , h — 1) = rank A(i + 1, . . . , h — 1, h) = rank A(i, . . . , h) — 
rank A(i + 1, . . . , h — 1) + 1. Let Ii and 7j + i be the lexicographically minimal fc-subsets with respect to 
the total orders <i and <i+i, such that Aj i (A) ^ and A/ i+1 (A) ^ 0. Since rank A(i, i + 1, ... ,h — l) = 
vankA(i, . . . , h), we have h h. And since rank A(i + 1, . . . , h — 1, h) = rank A(i + 1, . . . , h — 1) + 1, we 
have h G h+i- We now claim that i G Ii. Otherwise, by the definition of Grassmann necklace, J^+i = Ii, 
which contradicts the fact that rank A(i, i + 1, . . . , h — 1) = rank A(i + 1, . . . , h — 1) + 1. Therefore the 
claim holds, and by Definition 2.14, we must have ir(h) = i. □ 
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9. SOLITON GRAPHS AND GENERALIZED PLABIC GRAPHS 

The following notion of soliton graph forgets the metric data of the contour plot, but preserves the 
data of how line-solitons interact and which exponentials dominate. 

Definition 9.1. Let A G Grk >n and consider a generic contour plot Ct(uA) for some time t. Color a 
trivalent vertex black (respectively, white) if it has a unique edge extending downwards (respectively, 
upwards) from it. We preserve the labeling of regions and edges that was used in the contour plot: we 
label a region by Ej if the dominant exponential in that region is AjEj, and label each line-soliton 
by its type (see Lemma 7.5). We also preserve the topology of the graph, but forget the metric 
structure. We call this labeled graph with bicolored vertices the soliton graph Gt (uAj- 

Example 9.2. We continue Example 7.9. Figure 4 contains the same contour plot Ct{uA) as that at the 
left of Figure 3. One may use Lemma 7.5 to label all regions and edges in the soliton graph. After 
computing the Pliicker coordinates, one can identify the singular solitons, which are indicated by the 
dotted lines in the soliton graph. 




FIGURE 4. Example of a contour plot C t (uA), its soliton graph C = G t (uA), and its 
generalized plabic graph Pl{C). The parameters used are those from Example 7.9. In 
particular, [k\, . . . , Kg) — (—5, —3, —2, —1, 0, 1, 2, 3, 4), and n — (6, 7, 1, 8, 2, 3, 9, 4, 5). 

We now describe how to pass from a soliton graph to a generalized plabic graph. 

Definition 9.3. A generalized plabic graph is an undirected graph G drawn inside a disk with n boundary 
vertices labeled {1, . . . , n}. We require that each boundary vertex i is either isolated (in which case it is 
colored with color 1 or —1), or is incident to a single edge; and each internal vertex is colored black or 
white. Edges are allowed to cross each other in an X-crossing (which is not considered to be a vertex). 

By Theorem 8.1, the following construction is well-defined. 

Definition 9.4. Fix a positroid stratum S^- of Grk, n where 7r = (tt,coI). To each soliton graph C 
coming from a point of that stratum we associate a generalized plabic graph Pl{C) by: 

• embedding C into a disk, so that each unbounded line-soliton of C ends at a boundary vertex; 

• labeling the boundary vertex incident to the edge with labels i and 7r(«) by 7r(z); 

• adding an isolated boundary vertex labeled h with color 1 (respectively, —1) whenever h G J 
for each region label Ej (respectively, whenever h ^ J for any region label Ej); 

• forgetting the labels of all edges and regions. 

See Figure 4 for a soliton graph C together with the corresponding generalized plabic graph Pl(C). 
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Definition 9.5. Given a generalized plabic graph G, the trip Ti is the directed path which starts at the 
boundary vertex i, and follows the "rules of the road" : it turns right at a black vertex, left at a white 
vertex, and goes straight through the A-crossings. Note that Ti will also end at a boundary vertex. If 
i is an isolated vertex, then T, starts and ends at i. Define no(i) = j whenever Ti ends at j. It is not 
hard to show that ttq is a permutation, which we call the trip permutation. 

We use the trips to label the edges and regions of each generalized plabic graph. 

Definition 9.6. Given a generalized plabic graph G, start at each non-isolated boundary vertex i and 
label every edge along trip Ti with i. Such a trip divides the disk containing G into two parts: the part 
to the left of Ti, and the part to the right. Place an i in every region which is to the left of Ti. If h is 
an isolated boundary vertex with color 1, put an h in every region of G. After repeating this procedure 
for each boundary vertex, each edge will be labeled by up to two numbers (between 1 and n), and each 
region will be labeled by a collection of numbers. Two regions separated by an edge labeled by both i 
and j will have region labels S and (S \ {i}) U {j}. When an edge is labeled by two numbers i < j, we 
write on that edge, or {i, j} or {j, i} if we do not wish to specify the order of i and j. 

Although the following result was proved for irreducible cells of (Grk, n )>o, the same proof holds for 
arbitrary positroid strata of Gru,n- 

Theorem 9.7. [20, Theorem 7.6] Consider a soliton graph C — G t {uA) coming from a point A of a 
positroid stratum S v ., where 7r = (n,col). Then the trip permutation of Pl(C) is ir, and by labeling 
edges of Pl(C) according to Definition 9.6, we will recover the original edge and region labels in C. 

We invite the reader to verify Theorem 9.7 for the graphs in Figure 4. 

Remark 9.8. By Theorem 9.7, we can identify each soliton graph C with its generalized plabic graph 
Pl(C). From now on, we will often ignore the labels of edges and regions of a soliton graph, and simply 
record the labels on boundary vertices. 



Consider a matroid stratum Sm contained in the Deodhar component Sd, where D is the corre- 
sponding or Go-diagram. From Definition 7.1 it is clear that the contour plot associated to any A e Sm 
depends only on Ai, not on A. In fact for t < a stronger statement is true - the contour plot for 
any A e Sm C Sd depends only on D, and not on A4. In this section we will explain how to use D to 
construct first a generalized plabic graph G-(D), and then the contour plot C t {M) for i<0. 

10.1. Definition of the contour plot for (<0. Recall from (7.1) the definition of fM(x,y,t). To 
understand how it behaves for t -C 0, let us rescale everything by t. Define x = f and y = | , and set 



that is, KiX + n 2 y + n^t = t(f>i(x, y). Note that because t is negative, x and y have the opposite signs of 
x and y. This leads to the following definition of the contour plot for t <C 0. 

Definition 10.1. We define the contour plot C_ 00 (A / J) to be the locus in R 2 where 



Remark 10.2. After a 180° rotation, C_ co (A'i) is the limit of C t (uA) as t — >• — oo, for any A e Sm- 
Note that the rotation is required because the positive x-axis (respectively, y-axis) corresponds to the 
negative x-axis (respectively, y-axis). 



10. The contour plot for K<0 



<t>i{x,V) = KiX + n 2 y + k\ , 
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Definition 10.3. Define Vi^, m to be the point in M 2 where (f>i(x,y) = <fre(x,y) — <j> m (x,y). A simple 
calculation yields that the point Vi^_ m has the following coordinates in the xy-plane: 

Some of the points v it £ m G R 2 correspond to trivalent vertices in the contour plots we construct; 
such a point is the location of the resonant interaction of three line-solitons of types [i,£], [i,m] and 
[i, m] (see Theorem 10.6 below). Because of our assumption on the genericity of the K-parameters, those 
points are all distinct. 

10.2. Main results on the contour plot for t<0. The results of this section generalize those of 
[19, Section 8] to a soliton solution coming from an arbitrary point of the real Grassmannian (not just 
the non-negative part). We start by giving an algorithm to construct a generalized plabic graph G_(D), 
which will be used to construct C_ oc (A / J). Figure 5 illustrates the steps of Algorithm 10.4, starting from 
the Go-diagram of the Deodhar component Sd where D is as in the upper left corner of Figure 5. 

Algorithm 10.4. From a Go-diagram D to G-(D): 

(1) Start with a Go-diagram D contained in a k x (n — k) rectangle, and replace each O, 9, and 
blank box by a cross, a cross, and a pair of elbows, respectively. Label the n edges along the 
southeast border of the Young diagram by the numbers 1 to n, from northeast to southwest. 
The configuration of crosses and elbows forms n "pipes" which travel from the southeast border 
to the northwest border; label the endpoint of each pipe by the label of its starting point. 

(2) Add a pair of black and white vertices to each pair of elbows, and connect them by an edge, as 
shown in the upper right of Figure 5. Forget the labels of the southeast border. If there is an 
endpoint of a pipe on the east or south border whose pipe starts by going straight, then erase 
the straight portion preceding the first elbow. If there is a horizontal (respectively, vertical) 
pipe starting at i with no elbows, then erase it, and add an isolated boundary vertex labeled i 
with color 1 (respectively, —1). 

(3) Forget any degree 2 vertices, and forget any edges of the graph which end at the southeast 
border of the diagram. Denote the resulting graph G-(D). 

(4) After embedding the graph in a disk with n boundary vertices (including isolated vertices) we 
obtain a generalized plabic graph, which we also denote G_(D). If desired, stretch and rotate 
G_(D) so that the boundary vertices at the west side of the diagram are at the north instead. 

Remark 10.5. If there are no black stones in D, then this algorithm reduces to [20, Algorithm 8.7]. 
In this case, by [20, Theorem 11.15], the Pliickcr coordinates corresponding to the regions of G-(D) 
include the set of minors J described in Theorem 5.12. In particular, the set of Plucker coordinates 
labeling the regions of G_ (D) comprise a positivity test for Sd ■ 

The following is the main result of this section. 

Theorem 10.6. Choose a matroid stratum Sm an d let Sd be the Deodhar component containing Sm- 
Recall the definition of n(D) from Definition 4.10. Use Algorithm 10.4 to obtain G-(D). Then G-(D) 
has trip permutation n(D), and we can use it to explicitly construct C- tx (A4) as follows. Label the 
edges of G-(D) according to the rules of the road. Label by 

Vi,£,m each trivalent vertex which is 
incident to edges labeled [i,£], [i,m], and [£, m], and give that vertex the coordinates (x,y) — (k^k^ + 
^iK-m "I" K £ K mj — («i + k £ + K m )). Replace each edge labeled which ends at a boundary vertex by an 
unbounded line-soliton with slope Ki + Kj. (Each edge labeled [i,j] between two trivalent vertices will 
automatically have slope Ki + Kj.) In particular, C_ 00 (A / J) is determined by D. Recall from Remark 
10.2 that after a 180° rotation, C_ co (A / J) is the limit of C t (uA) as t — > — oo, for any A e Sm- 

Remark 10.7. Since the contour plot C_ 00 (A / I) depends only on D, we also refer to it as C-oo(D). 
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Figure 5. Construction of the generalized plabic graph G-(D) associated to the Go- 
diagram D. The labels of the regions of the graph indicate the index sets of the 
corresponding Pliicker coordinates. Using the notation of Definition 4.10, we have 
tt(D) = vw- 1 = (5, 7, 1, 6, 8, 3, 4, 2). 



Remark 10.8. The results of this section may be extended to the case t ^> by duality considerations 
(similar to the way in which our previous paper [20] described contour plots for both t <C and t 3> 0). 
Note that the Deodhar decomposition of Grk,n depends on a choice of ordered basis (ei, . . . , e n ). Using 
the ordered basis (e„, . . . , e\) instead and the corresponding Deodhar decomposition, one may explicitly 
describe contour plots at t » 0. 

Remark 10.9. Depending on the choice of the parameters Kj, the contour plot C-oo(D) may have a 
slightly different topological structure than the soliton graph G_(D). While the incidences of line- 
solitons with trivalent vertices are determined by G_(D), the locations of X-crossings may vary based 
on the Kj's. More specifically, changing the «j's may change the contour plot via a sequence of slides, 
see Section 11. 

Our proof of Theorem 10.6 is similar to the proof of [19, Theorem 8.9]. The main strategy is to use 
induction on the number of rows in the Go-diagram D. More specifically, let D' denote the Go-diagram 
D with its top row removed. In Lemma 10.11 we will explain that G-(D') can be seen as a labeled 
subgraph of G-(D). In Theorem 10.14, we will explain that there is a polyhedral subset of C_oo(-D) 
which coincides with C- oa {D'). And moreover, every vertex of C_ 00 (D') appears as a vertex of C_ oc (£'). 
By induction we can assume that Theorem 10.6 correctly computes C_ oc (Z3'), which in turn provides 
us with a description of "most" of C_oo(D), including all line-solitons and vertices whose indices do 
not include 1. On the other hand, Theorem 8.1 gives a complete description of the unbounded solitons 
of both C-oo(D') and C-oo(D) in terms of ft(D') and tt(D). In particular, C-oo(D) contains one more 
unbounded soliton at y than does C_ 00 (Z?). This information together with the resonance property 
allows us to complete the description of C_ oc (Z3) and match it up with the combinatorics of G-(D). 
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Lemma 10.10. The generalized plabic graph G-(D) from Algorithm 10.4 has trip permutation tt(D). 

Proof. If we follow the rules of the road starting from a boundary vertex of G_ (D) , we will first follow 
a "pipe" southeast (compare the lower left and the top middle pictures in Figure 5) and then travel 
straight west along the row or north along the column where that pipe ended. Recall from Definition 
4.10 that 7t(.D) = vw~ . Noting that we can read off v and w from the pipes in the top middle picture 
of Figure 5, we see that following the rules of the road has the same effect as computing vw~ 1 . □ 

The next lemma explains the relationship between G_(D) and G-(D'), where D' is the Go-diagram 
D with the top row removed. It should be clear after examining Figure 6. 





o 


o 




• 








• 


o 




4 3 




o 


o 


5 


8 


7 


6 




• 








• 


o 




3 

4 




o 


o 


5 




[2,8] [4,7] [3,6] [1,3] 



7 6 



5 4 3 2 





[5,8] [4,7] [3,6] [2,3] 



[2,7] - 

[4,6] 

[5,8] 




• 


O 






O 


o 


8 7 


6 




o 


o 


8 7 


6 



5 4 



[4,6] - 
[5,8] - 



[5,8] [4,6] 



Figure 6. Inductive construction of the generalized plabic graph G-(D) associated 
to the Go-diagram D, cf. Figure 5. 



Lemma 10.11. Let D be a Go-diagram with k rows and n—k columns, and let G = G- (D) be the edge- 
labeled plabic graph constructed by Algorithm 10.4. Form a new Go-diagram D' from D by removing 
the top row of D; suppose that I is the sum of the number of rows and columns in D' . Let G' be the 
edge-labeled plabic graph associated to D' , but instead of using the labels {1, 2, . . . ,£}, use the labels 
{n — I + 1, n — £ + 2, . . . , n}. Let h denote the label of the top row of D. Then G' is obtained from G 
by removing the trip T/j starting at h and all edges to its right which have a trivalent vertex on Th- 

From now on, we will assume without loss of generality that i\ = 1 is a pivot for A € Sd- 

Definition 10.12. Let A4 be a matroid on [n] such that 1 is contained in at least one base. Let M! be 
the matroid { J \ {1} | 1 € J and J € M}. 

Using arguments similar to those in the proof of Theorem 5.6, one can verify the following. 
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Contour plot for \\ 



C_oo(M') 




Figure 7. The contour plot C-^M') within the contour plot C_ 00 (.M). 



Lemma 10.13. If A G Sm C Sd is in row-echelon form and A! is the span of rows 2, 3, . . . , fc in 
A G SVi C Grk t n, then A' G Sjw C Sd>, where D' is obtained from D by removing its top row. 

The following result is a combination of [19, Theorem 8.17] and [19, Corollary 8.18]. Although in 
[19] the context was A G (Grk, n )>o and in this paper we are allowing A G GVfc in , the proofs from [19] 
hold without any modification. See Figure 7 for an illustration of the theorem. 

Theorem 10.14. [19] Let M. be a matroid such that 1 is contained in at least one base. Then there is 
an unbounded polyhedral subset TZ of C-oo(M) whose boundary is formed by line-solitons, such that 
every region in TZ is labeled by a dominant exponential Ej such that 1 G J. In TZ, C- 00 (Ai) coincides 
with C~oo(A4'). Moreover, every region of C-^M.') which is incident to a trivalent vertex and labeled 
by Eji corresponds to a region of C_ 00 (A / 1) which is labeled by T£/'u{i}- 

In particular, the set of trivalent vertices in C-oo(M) is equal to the set of trivalent vertices in 
C-oo(M) together with some vertices of the form V\ & c . These vertices comprise the vertices along the 
trip Ti (the set of line-solitons labeled [l,j] for any j). In particular, every line-soliton in C_ oc (A / J) 
which was not present in C—oo(M.') and is not on T\ must be unbounded. And every new bounded 
line-soliton in C-oo(-M) that did not come from a line-soliton in C_ 00 (A / I / ) is of type [1, j] for some j. 

We now prove Theorem 10.6, using the characterization of unbounded line-solitons in Theorem 8.1. 

Proof. Choose A in the Deodhar component Sd- Let Ai be the matroid such that A G Sm- We 
will prove Theorem 10.6 using induction on the number of rows of A. Using the notation of Definition 
10.12 and Lemma 10.13, we have that A' G Sm' C Sd>- 

By Theorem 10.14, the contour plot C_ 00 (.M) is equal to the contour plot C-oo(M') together with 
some trivalent vertices of the form i>i,6, c , all edges along the trip Ti, and some new unbounded line- 
solitons (which are all to the right of the trip Ti). By the inductive hypothesis, C^ ao {M') is constructed 
by Theorem 10.6; in particular, Algorithm 10.4 produces a (generalized) plabic graph which describes 
the trivalent vertices of C_ 00 (A / 1') and the interactions of all line-solitons at trivalent vertices. 

Using Lemma 10.10 and Theorem 8.1, we see that Algorithm 10.4 produces a (generalized) plabic 
graph whose labels on unbounded edges agree with the labels of the unbounded line-solitons for the 
contour plot C- OD (M) of any A G Sd- The same is true for A' G Sd 1 - 

By Lemma 10.11, the plabic graph G which Algorithm 10.4 associates to D is equal to G' together 
with the trip Ti starting at 1 at some new line-solitons emanating right from trivalent vertices of T\ . 
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We now characterize the new vertices and line-solitons which C-oo{M.) contains, but which C-oo(M') 
did not. We claim that the set of new vertices is precisely the set of Vi & c (where 1 < b < c), such that 
either c — > b is a nonexcedance of 7r = ir(A4), or c — >• 6 is a nonexcedance of 7r' = 7r(.M'), but not both. 
Moreover, if c — > & is a nonexcedance of 7r, then i>i.f, iC is white, while if c — > 6 is a nonexcedance of 7r', 
then vi^.c is black. The proof is identical to that of the same claim in the proof of [20, Theorem 8.8]. 

Now, if one analyzes the steps of Algorithm 10.4 (see in particular the second and third diagrams 
in Figure 5), it becomes apparent that the above description also characterizes the set of new vertices 
which the algorithm associates to the top row of the Go-diagram D. In particular, the nonexcedances 
of the corresponding permutation 7r correspond to the vertical edges at the top of the second and third 
diagrams; when one labels these edges using the rules of the road, each edge gets the label [b, c], where 
b comes from the label of its pipe, and c comes from the label of its column (shown at the bottom of the 
second diagram). The nonexcedances of n' are labeled in the same way but come from vertical edges 
which are present in the second row of D. Therefore each new trivalent vertex in the top row gets the 
label ui,6 )C where b and c are as above, and where c — ► b is a nonexcedance of precisely one of 7r and it' . 

Finally, we discuss the order in which the vertices v± t b,c occur along the trip T\ in the contour plot. 
First note that the trip T\ starts at y <C and along each line-soliton it always heads up (towards 
y ^> 0). This follows from the resonance condition (see e.g. [20, Figure 9] and take i — 1). Therefore 
the order in which we encounter the vertices V\ & iC along the trip is given by the total order on the 
y-coordinates of the vertices, namely K\ + Kb + K c . 

We now claim that this total order is identical to the total order on the positive integers 1 + b + c 
that is, it does not depend on the choice of k^'s, as long as K\ < • • • < K n . If we can show this, then 
we will be done, because this is precisely the order in which the new vertices occur along the trip T\ in 
the graph G_(L). 

To prove the claim, it is enough to show that among the set of new vertices v\^,c, there are not two 
of the form v\^j and v\ j.k where i < j < k < I. To see this, recall that the indices b and c of the new 
vertices i>i.f, iC can be read off from the second and third diagrams illustrating Algorithm 10.4: c will 
come from the bottom label of the corresponding column, while b will come from the label of the pipe 
that vi.b^c lies on. Therefore, if there are two new vertices v\^j and vij,fc, then they must come from 
a pair of pipes which have crossed each other an odd number of times, as in Figure 8. 



Vi,u 

~> Vljjc- 



k' k 



] J 
Vu,i k J Vijf 



-(H)- 



k' k 



Figure 8 



Note that the second diagram of Figure 5 depicts a "pipe dream" (or "wiring diagram" ) encoding 
the distinguished subexpression v of a reduced expression w. If two pipes pass over each other in a 
given box we will say that they cross at that box, while if two pipes pass through the same box without 
crossing, we will say that they kiss at that box. Let us now follow a pair of pipes from southeast to 
northwest. The property of v being distinguished mean that two pipes starting at i and j must not 
kiss each other after having crossed each other an odd number of times. 
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Assume that Algorithm 10.4 produces two vertices i>i,i/ and t>i,j,fc where i < j < k < £. Choose such 
a pair of vertices which minimizes \£ — k\. We consider two cases, based on whether V\j f. is black or 
white. In the first case (see the left of Figure 8), since Vij t k is black, its pipe j will continue west from 
i>i,j,fc and must eventually turn up, at some column k' such that k < k' < £■ But then Algorithm 10.4 
produces another vertex i>i,j,fe' such that i < j < k' < t, so this vertex together with v\^j form a pair 
of vertices where \t — k'\ < \l — k\, contradicting our assumption of minimality of \£ — k\. 

In the second case (see the right of Figure 8), since Vx j,k is white, there is another black vertex f i j'.fc 
to its left in the same box b, whose pipe starts at j'. Because v is distinguished, j' must be greater 
than j. (Otherwise the pipes starting at j and j' would cross each other an odd number of times and 
then kiss at box b.) Now since is black, its pipe must travel west from it and eventually turn up, 

at some column k' such that f < k' < I. But then Algorithm 10.4 produces another vertex Viji^i such 
that i < j' < k' < £. But now we have a pair of vertices Uj. » £ and V\ji^i such that i < j' < kl < £ 
where \£ — k'\ < \£ — k\. This contradicts our assumption of minimality of \£ — k\, and completes the 
proof of the claim. 

Finally, using Definition 10.3 for the vertex Uj^ )TO , we obtain the contour plot from G by giving the 
trivalent vertices the explicit coordinates from Theorem 10.6. □ 

11. X-CROSSINGS, SLIDES, AND CONTOUR PLOTS 

In this section we discuss how our choice of the parameters Ki may affect the topology of the contour 
plot C-oo(D) (and hence Ct(uA) for t and A € Sd), namely, by changing the locations of the X- 
crossings. See Remark 10.9. We also discuss the relation between X-crossings and Pliicker coordinates. 

11.1. Slides and the topology of contour plots. The following definition will be useful for under- 
standing the dependence of the contour plot on the k^'s. 






o-o- 
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FIGURE 9. Some slides involving white X-crossings. These contour plots correspond 
to the same Le-diagram D with tt(D) — (5, 3, 2, 1, 4), but they differ from G-{D). 



Definition 11.1. Consider a generalized plabic graph G with at least one A-crossing. Let i; a ,b.c be a 
trivalent vertex (with edges labeled [a, b], [a, c], and [&, c]) which has a small neighborhood N containing 
one or two X-crossings with a line labeled but no other trivalent vertices or X-crossings. Here 

{a, b, c} and {i, j} must be disjoint. Then a slide is a local deformation of the graph G which moves the 
line [i,j] so that it intersects a different set of edges of v a b >c , creating or destroying at most one region 
in the process. 

See Figure 9 for examples. Recall the notions of black and white X-crossings from Definition 7.7. 

Remark 11.2. Theorem 10.6 determines everything about the combinatorics and topology of the contour 
plot C_oo(D) except for which pairs of line-solitons form an X-crossing. Therefore if one deforms the 
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Figure 10. Contour plots C t {uA) constructed using the same t and A £ Sd C Gr^g 
but with different choices of the K-parameters. The left plot uses (Ki,...,/tg) = 
(-3.5, -2,-1, 0, 0.5, 1, 2, 5) while the right one uses (-3.5, -2,-1, 0, 0.5, 1, 2.5, 3). This 
affects the location of the [4, 7] line-soliton. In the middle we have the generalized plabic 
graph G_ (D) using the Go-diagram D of Figure 5. 

parameters k^, the only way that the contour plot can change so as to change the topology is via a 
sequence of slides. 

See Figure 10 for an example of two different contour plots associated to the same Go-diagram and 
element A G Gr^.%, but obtained using different choices of the K-parameters. The two contour plots 
differ by precisely one slide. For another example, compare Figure 4 to Figure 11. Both of them are 
based on the Go-diagram from Example 7.9 and the same matrix A. The only difference is the value 
of K\. Note that this affects the X-crossings formed by the unbounded [1,6] line-soliton, and that one 
contour plot can be obtained from the other via a sequence of three slides. 




Figure 11. A contour plot C t (uyt), soliton graph C — Gt{uA) and generalized plabic 
graph G-(D) coming from a Go-diagram where A E Sd- The ^-parameters are the 
same as those used for Figure 4 except that ki = —3.1 now, i.e. (ki, . . . ,kq) = 
(-3.1,-3,-2,-1,0,1,2,3,4). 

We now show that a slide on a contour plot preserves the number of black X-crossings. 

Theorem 11.3. Consider two contour plots C and C' (for the same A 6 Gr^^n and time t but for different 
K-parameters) which differ by a slide. Then C and C' have the same number of black X-crossings. 
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Proof. Suppose that C and C differ by a slide involving the trivalent vertex v a ^, c and the line-soliton 
for a < b < c and i < j, where the sets {a, b, c} and are disjoint. We assume that v a b, c is 

white. (The case where it is black is analogous.) There are five cases to consider: 
Case 1. i<a<j<b<c, which implies that m + Kj < n a + Kb < k c < K b . 
Case 2. i<a<b<j<c, which implies that (a.) Ki + Kj < K a + K b < n a + k c < Kb + K c , or 

(b.) K a + K b < Ki + Kj < K a + K c < K b + K c . 

Case 3. a<i<b<j<c, which implies that (a.) K a + Kb < Ki + Kj < K a + k c < Kb + k c , or 

(b.) K a + K b < K a + K c < Ki + K-j < K b + K c - 

Case 4. a<i<b<c<j, which implies that (a.) K a + K b < K a + k c < Ki + Kj < K b + k c , or 

(b.) K a + K b < K a + K c < K b + K c < Ki + Kj. 

Case 5. a<b<i<c<j, which implies that K a + K b < K a + k c < K b + k c < Ki + Kj. 
(Note that any other ordering on a,b,c,i,j, such &si<j<a<b<c, would imply that there are no 
black X-crossings involving the edges incident to v a ,b,c and the [i, j] soliton.) 




Figure 12. Various types of X-crossings involving the line-solitons incident to v a ,b,c 
and the [i,j] line-soliton. The top row shows Cases 1, 2a, 2b, and 3a from left to right, 
while the bottom row shows Cases 3b, 4a, 4b, and 5 from left to right. 

Consider Case 1. Recall that "slope" of the line-soliton - that is, the tangent of the angle 
measured counterclockwise from the positive y-axis to the [i, j] line-soliton - is equal to Ki+Kj . Therefore 
from the order on the slopes, the soliton may intersect either the [a, c] soliton or both the [a, b] 
and [b, c] solitons, as in the top-left diagram of Figure 12. The black X-crossings are denoted by a solid 
black square. In both cases, precisely one of the intersections is a black X-crossing. The other cases 
are similar - see Figure 12. □ 

Remark 11.4. In fact one can show that the slides from Cases 3a and 3b in Figure 12 are impossible at 
( < 0. More specifically, it is impossible for the [i,j] line-soliton to intersect the [b, c] line-soliton. To 
show this, one may compute the coordinates (x v , y v ) of the trivalent vertex v where the [a, b], [a, c], and 
[b, c] solitons intersect. Then one can show that the intersection of the [i,j] soliton and the line y = y v 
has ^-coordinate which is strictly less than x v . 
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11.2. Slides and Plucker coordinates. In [20, Theorem 9.1], we proved that the presence of X- 
crossings in contour plots at \t\ » implies that there is a two-term Plucker relation. 

Theorem 11.5. [20, Theorem 9.1] Suppose that there is an X-crossing in a contour plot Ct(uA) for some 
A e Grk.n where \t\ ^S> 0. Let h, I2, h, and I4 be the /c-element subsets of {1, ... , n} corresponding to 
the dominant exponentials incident to the X-crossing listed in circular order. 

• If the X-crossing is white, we have A Il (A)A l3 (A) = A l2 (A)A h (A). 

• If the X-crossing is black, we have A Il (A)A h (A) = -A l2 (A)A Ii (A). 

The following corollary is immediate. 

Corollary 11.6. If there is a black X-crossing in a contour plot at i < or t > 0, then among the 
Plucker coordinates associated to the dominant exponentials incident to that black X-crossing, three 
must be positive and one negative, or vice-versa. 

Corollary 11.7. Let D be a J-diagram, that is, a Go-diagram with no black stones. Let A e Sd and 
(<0. Choose any m < • • • < n n . Then the contour plot C((ua) can have only white X-crossings. 

Proof. From Theorem 10.6, it follows that the contour plot C-^ua) has no dependence on the signs 
of the Plucker coordinates of A. (In fact it has no dependence on A, only on the Deodhar stratum 
Sd containing A.) Since I? is a J-diagram, we can choose an element A' £ Sd H (Gr^„)> , and 
C-oo(ma) = C-co("4')- But now since the Plucker coordinates of A' are all non-negative, by Theorem 
11.5, there cannot be any black X-crossings in the contour plot. □ 

Lemma 11.8. Consider two contour plots for A e Grk.n which differ by a single slide. Let J and 
J' denote the two sets of Plucker coordinates corresponding to the dominant exponentials in the two 
contour plots. Then from the values of the Plucker coordinates in J , one can reconstruct the values of 
the Plucker coordinates in J' , and vice- versa. 

Proof. By Theorem 11.5, the four Plucker coordinates incident to an X-crossing satisfy a "two-term" 
Plucker relation. Now it is easy to verify the lemma by inspection, since each slide only creates or 
removes one region, and there is a dependence among the Plucker coordinates labeling the dominant 
exponentials. The reader may wish to check this by looking at the first and second, or the second and 
third, or the third and fourth contour plots in Figure 9. □ 

Corollary 11.9. Let D be a J-diagram, such that Sd C Grk.n- Let C-oo(D) and C / _ ca (D) be two 
contour plots defined using two different sets of parameters K\ < ■ ■ ■ < K n and k\ < ■ ■ ■ < n' n . Let J 
and J' be the fc-element subsets corresponding to the dominant exponentials in C_ co (D) and C_ !X (D). 
If Aj(A) > for each I & J, then Ai{A) > for each I e J'. In particular, if J' is a positivity test 
for Sd then so is J'. 

Proof. One may use a continuous deformation of the parameters to get from n x < ■ ■ ■ < n n to 
one deforms the parameters the contour plot will change by a sequence of slides. At 
each step along the way, the contour plot will contain only white X-crossings (by Corollary 11.7). By 
Lemma 11.8, if we know the values of the Plucker coordinates labeling dominant exponentials before 
a slide, then we can compute the Plucker coordinates labeling dominant exponentials after a slide. 
Moreoever, since this computation involved only two-term Plucker relations and all the X-crossings are 
white, the positivity of the Plucker coordinates in J implies the positivity of the Plucker coordinates 
in J'. □ 

12. The regularity problem for KP solitons 

In this section, we first discuss the regularity of KP solitons. Given a soliton solution ua coming 
from an element A e Gr k ^ n , we show that if Ua{x,u, t) is regular for f« 0, then in fact A must lie 
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in the totally non-negative part (Grk, n )>o of the Grassmannian. We then discuss the uniqueness (and 
lack thereof) of the pattern when the soliton solution is not regular. 
Our main theorem is the following. 

Theorem 12.1. Fix parameters Ki < ■ ■ ■ < K n and an element A £ Gr kn . Consider the corresponding 
soliton solution ua(x,U, t) of the KP equation. This solution is regular at t <C if and only if A £ 
(Gr k>n )> . Therefore this solution is regular for all times t if and only if A £ (Gr k ^ n )> . 

We will prove Theorem 12.1 in Section 12.2, after establishing some results on black X-crossings. 

12.1. Lemmas on black X-crossings. Recall from Section 10.1 that (f>i(x, y) — nix + nfy + nf. The 
following lemma is easy to check. 

Lemma 12.2. For 1 < i < j < n, let be the line in the xy -plane where <pi(x,y) = (f>j(x,y). For 
i < j < k < £, let bij t k,t be the point where the lines Li k and Ljg intersect. Then Lij has the equation 

x + («j + Kj)y + (k? + KiKj + Kj) = 0, 
and the point bij. k .g = (b* j k g, b\ • k £ ) has the coordinates 

2 22 222 22 

h,j,k,e = — T Z 

K( Kj -\- K k Kg 

K? KiK k "i" Kj "i" KjKg fc'fc ~t" Kg 
' Ki Kj ~\~ Kk Kg 

Lemma 12.3. Consider the point bij^,e where 1 ^ {i,j,k,£}. Then at this point we have <j)\ < fa — fa 
and fa < <pj = <f>£. 

Proof. By definition of h,j t k,e we have that at this point fa = fa and (f>j = fa. So we just need to 
show that at bij tkt g, fa < fa and fa < fa:. A calculation shows that 4>i{bi,j,k,i) — fa{bij.k,i) is equal to 

(K k - Ki)(Kj - Kl)[(/Cj - /Ci)(/Cj - Kj+Kj- K k ) + (Kj - Kjjijig - K k )] 

Kj Ki -f- Kg K k 

and (pj(bij tk ,e) ~ fa(h,j,k,e) is equal to 

(Kj - Kl)(Kj - Kl)[(Kj - Kl)(K 3 - Kj + Kj - K k ) + (k £ - K k )(K k - Kj)] 

Kj — Ki + Kg — K k 

Because K\ < Ki < Kj < K k < Kg, we can readily verify that the above quantities are positive. □ 

Remark 12A. Lemma 12.3 will be instrumental in proving Proposition 12.5 below regarding black X- 
crossings. Note that if in the lemma we took the order i < k < j < £ or i < j < £ < k then our proof 
would not work. So Proposition 12.5 does not necessarily hold for white X-crossings. 

Proposition 12.5. Use the hypotheses and notation of Theorem 10.14. Then every black X-crossing of 
C-oo(-M') remains a black X-crossing in C-oo(M); and each region in C_ 00 (A1 / ) which is incident to a 
black X-crossing and is labeled by Ej> corresponds to a region of C_ oc> (M) which is labeled by -Ej'u{i}- 

Proof. Consider a black X-crossing b a ^, c ,d of C-^M.') in which the linc-solitons [a,c] and [6, d] 
intersect (here a < b < c < d). Since this is taking place in C_ 00 (A / J'), 1 ^ {a, b, c, d}. The four regions 
Ri,R2,Rz,Ra incident to b a ^, Ci d are labeled by Ej 1 , Ej 2 , Ej 3 , Ej 4 . In particular, this means that at 
region R\, J\ is the subset {ji, . . . , j k -i} of M' which minimizes the value 6j 1 +■ ■ • + dj k _ 1 - Without loss 
of generality we can assume that a £ J\. But then by Lemma 12.3, there is a neighborhood TV of b a .b,c.d 
where <f>\ is less than <f) a . It follows that in N fl Ri, J\ U {j k = 1} is the subset of M that minimizes 
the value 6j 1 + ■ ■ ■ + 9j k . Therefore the region Ri of C_oo (M r ) which is labeled by E Jl corresponds to a 
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region of C- ao (AA) which is labeled by £'j 1 u{i}- Similarly for R2, R3, and R4. In particular, the black 
X-crossing from C_ 00 (A / {') will remain a black X-crossing in C-oo(M). □ 
Recall the notion of a slide from Definition 11.1. 

Proposition 12.6. Choose a Go-diagram D such that Sd C Grk >n - Let K\ < ■ ■ ■ < n n and n[ < ■ ■ ■ < n' n 
be two choices of parameters, and let C_oo(D) and C_ oa (D) be the corresponding contour plots. Then 
if C-^D) has r black X-crossings, then C'_ QO (D) has r black X-crossings. 

Proof. By Remark 11.2, the two contour plots differ by a series of slides. And by Theorem 11.3, each 
slide preserves the number of black X-crossings. □ 

Theorem 12.7. If D is a Go-diagram with at least one black stone, then the contour plot C-oo(D) 
contains a black X-crossing. 

Proof. Let i denote the bottom-most row of D which contains a black stone. Choose A £ Sd and 
put it in row-echelon form; let A' denote the span of rows i, i + 1, . . . , k of A. So A' G Sd> C Gr fc _ i+1 ,„, 
where D' is the Go-diagram obtained from rows i, i + 1, . . . , k of D. Then by Proposition 12.5, if we 
can show that the contour plot C-^D') contains a black X-crossing, then C_ 00 (Z?) must also contain 
a black X-crossing. 

Our goal now is to show that there is a choice of the /{-parameters such that C_ 00 (D / ) contains a 
black X-crossing. If we can show this, then by Proposition 12.6, we will be done. 
Note that for t = — 1, we have the following. 

(i) If i < j < k, then the y-coordinate yij^k of the trivalent vertex Vij^ where the [j, k] and 

[i, k] solitons meet is: 



(ii) If i < j < k < i, then the y-coordinate yij,k,i of an X-crossing between the [i,k] and [j,£] 



Consider the left-most black stone b in D' . Let [i, b] and [a, j] with i < a < b < j be the pair of lines 
in G-(D') which cross at this black stone. Then there are two cases: 

(a) There is no empty box to the left of b in D' , and so there is an unbounded [i, £>]-soliton at y 3> 
in the corresponding contour plot. Because b is a black stone, the [i, 6]-soliton must have a 
trivalent vertex v^bji at one end, where j' > b. Additionally, [a,j] is an unbounded soliton at 
y <C 0, and it has a trivalent vertex Wi', a .j at one end, where i! < a. See Figure 13. 



Ui,j,k — K i T K j + K k- 
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Figure 13 
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If we can choose the K-parameters such that yi', a ,j > Vi,a,b,j > Vi,b,j' then it follows that 
there is an intersection of the [a, j] and [i, b] line-solitons in the contour plot. 
One simple choice is to require that 

(12.1) Kj = —Ki > and Kb = — n a > 0; and also 

(12.2) Ki' > — Ki and Kj' < —Kj. 

By (12.1), we have J/i, ,6,j = 0- By (12.2), together with K a > /t,-/ and Kb < Ky , we have that 

Vi',a,j = Ki' ~\~ K a -\- Kj > > Ki -\- Kij -\- Kj' = yi^bjf . 

One concrete choice of parameters satisfying the required inequalities is (k^, k,/, K ai k^/, k^) = 
(— 4r, — 2r, — r, r, 2r, 4r) where r > 0. 

(b) The second case is that there is an empty box to the left of b in D' , and so the [i, b] line-soliton 
has trivalent vertices at both ends. Figure 14 illustrates this situation. These vertices are the 
white vertex Vi^jn and the black vertex where i < i' < a < b < j' < j < j" . As before, 

[a,j] is an unbounded line-soliton at y <C which is incident to the trivalent vertex a j. Since 




Figure 14 
v i,b,j" is a white vertex, if we can show that 

yi,b,j" > Vi,a,b,j > Vi,b,j' \ and Vi'aJ > 2/i,a,b,j: 

then it follows that the line-solitons of type [<x,j] and [i, 6] intersect in the contour plot. 

As before, we choose the K-parameters so that (12.1) and (12.2) are satisfied. Then again we 
have yi, a ,b,j = 0, yi\ a ,j > 0, and yi,b,j' < 0. Note that any choice of Kj" > Kj gives yi,b,j" > 0, 
since Ki + Kb + Kj" > Ki + Kb + Kj = Kb > 0. 

This completes the proof. □ 

12.2. Positivity of dominant exponentials and the proof of Theorem 12.1. In this section we 
prove Theorem 12.8 below. Once we have proved it, the proof of Theorem 12.1 will follow easily. 

Theorem 12.8. Let A € Sd C Grk jn , where D is a J-diagram, and let t -C 0. If Aj(A) > for each 
dominant exponential Ej in the contour plot Ct(uA), then A £ (Grk, n )>o- In other words, the Pliicker 
coordinates corresponding to the dominant exponentials in Ct(uA) comprise a positivity test for Sd- 

Lemma 12.9. Theorem 12.8 holds for elements A € Gr l n . 
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Proof. Let A G Sd C Gri >n . If £> contains r empty boxes, then Sd has dimension r. Meanwhile, the 
element A will have precisely r+1 nonzero Pliicker coordinates. (We can normalize the lexicographically 
minimal one to be 1.) It is easy to see that G-{D) and hence C t {uA) will have r + 1 regions, each one 
labeled by a different dominant exponential corresponding to a Pliicker coordinate Aj(A) such that 
Aj(A) ^ 0. Therefore if each such Aj(A) > 0, then A e (Gn,„)> . □ 

Lemma 12.10. Let A G Sd C Grk, n - Then it is possible to choose ki < k 2 < • • • < k„ such that the 
unbounded linc-solitons at y <C in the corresponding contour plot C t {uA} (for any time t) appear in 
the same order as they do in the generalized plabic graph G_ (D) . 

Proof. Recall that in a contour plot, the unbounded line-solitons at y <C appear from left to 
right in increasing order of the slope Ki + Kj. While in G-(D), one may easily check that the unbounded 
line-solitons at y <C appear from left to right in increasing order of j. 

Now let us choose K\,. . . ,K n so that «j — Ki-i — r l for some constant r > 1. To prove the lemma, 
it suffices to prove that given two line-solitons [a, b] and [c, d] at y <C 0, where b < d, we have that 

(12.3) K a + Kb < k c + Kd, or cquivalently, — «b > n a — n c . 

Since a < b and c < d, we have a < d. By our choice of the Kj's, — Kb > r d . If a < c then 
K a - k c < 0, so (12.3) is obvious. On the other hand, if a > c, then « a — n c < r a + r a ~ 1 + • • • + 1 < r a+1 . 
And since a < d, equation (12.3) follows. □ 

We now prove Theorem 12.8. 

Proof. Our strategy is to use induction on the number of rows of A. Lemma 12.9 takes care of the 
base case of the induction. We suppose that A is in row-echelon form, and let A' be the element of 
Grk-i, n obtained from the bottom fc — 1 rows of A. Then A' G Sd' where D' is also a J-diagram (it is 
the restiction of D to its bottom k — 1 rows). Recall from Theorem 10.14 that "most" of the contour 
plot Ct(uA') is contained in the contour plot Ct(uA)- More precisely, every region of Ct{v,A>) which is 
incident to a trivalent vertex and labeled by Ej> corresponds to a region of C t (uA) which is labeled by 
Ej'u{i}- Because A is in row-echelon form with a pivot in row 1, Aj/ U { 1 j( j 4) = Aji(A'), so the fact 
that each A J , u ^ 1 y(A) > implies that Aj>(A') > 0. 

We now claim that all Pliicker coordinates corresponding to the dominant exponentials of the contour 
plot Ct(uA') are positive. To prove this, note that from C t {uA) we can in fact construct C t (uA')- all 
of the trivalent vertices of Ct(uA') are present in Ct(iiA), so it is just a matter of extending some line- 
solitons that were finite in Ct{uA) but arc unbounded in Ct(uA'). These line-solitons may create some 
new white X-crossings but cannot create black X-crossings, because D' is a J-diagram. If a single white 
X-crossing is created, then because three of its four regions are incident to a trivalent vertex, three 
of the four corresponding Pliicker coordinates are positive. But then by the two-term Pliicker relation 
relating the four Pliicker coordinates, the fourth Pliicker coordinate is positive as well. If multiple white 
X-crossings are created, then one can iterate the above argument, starting with a white X-crossing with 
three of its four regions incident to a trivalent vertex in the contour plot. This proves the claim. So by 
the inductive hypothesis, A' G (Gr k ^i^ n )> . 

Since A' G {Grk-i, n )>o, h follows that all the Pliicker coordinates labeling the regions of G-(D') are 
positive. And so all of the Pliicker coordinates labeling the regions of G_ (D) which correspond to the 
bottom fc — 1 rows of D are positive. (Recall again that Aj/ U { 1 j( j 4) = Aj>(A').) If we can show that 
the Pliicker coordinates labeling the regions of G_ (D) which come from the top row of D are positive, 
then by Remark 10.5, it will follow that A G (Grk, n )>o- 

By Lemma 12.10, we can deform the K-parameters so that the resulting contour plot C' t (uA) has its 
unbounded linc-solitons at y <C in the same order as those in G_ (D). Then the dominant exponentials 
at y <§C in C' t (ua) are precisely those of G_ (D) , which in turn come from the top row of D. By Corollary 
11.9, since the dominant exponentials of C t (u^) are positive, so are those of C' t (u a)- In particular, the 
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dominant exponentials of C' t (uA) at y <C are positive, so we can conclude that all of the Pliicker 
coordinates labeling the regions of G_(D) are positive. Therefore A € (Gr fcj „)> . □ 
Finally we are ready to prove Theorem 12.1. 

Proof. Recall the definition of UA(x,y,t) in terms of the r-function from Section 6.2. It is easy to 
verify that if A e (Grk. n )>o, then UA(x,y,t) is regular for all times t: the reason is that TA(x,y,t) is 
strictly positive, and hence its logarithm is well-defined. 

Conversely, let A e Grfc^, and suppose that UA{x,y,t) is regular for (<0. This means that the 
Pliicker coordinates Aj corresponding to the dominant exponentials in the contour plot C t {uA) must 
all have the same sign. Since the Grassmannian is a projective variety, we may assume that all of these 
Pliicker coordinates are positive. 

Let Sd be the Deodhar stratum containing A. If D has a black stone, then by Theorem 12.7, 
the contour plot C-oo(«a) contains a black X-crossing. But then by Corollary 11.6, two dominant 
exponentials incident to that black X-crossing must have opposite signs, which is a contradiction. 
Therefore we conclude that D has no black stones. It follows that the Deodhar diagram corresponding 
to I? is a J-diagram. But now by Theorem 12.8, it follows that A e {Grk t n)>o- 

Finally, note that if A e Gr^^ n and ua{x, y, t) is regular for all times t, then in particular it is regular 
for t <§; 0, so the arguments of the previous two paragraphs apply. Therefore A e (Grfc, n )> . □ 

Remark 12.11. Corollary 11.6 implies that there are singularities among the line-solitons forming a 
black X-crossing in a contour plot, and the singular solitons form a V-shape. 

Example 12.12. We revisit the example from Figures 5 and 10. Note that the contour plot at the left of 
Figure 10 is topologically identical to G_(D). The Go-diagram and labeled Go-diagram are as follows. 
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Recall from Theorem 5.6 that we associate a Pliicker coordinate A/ b to each box b of the Go-diagram, via 
lb = u m (w m ) _1 {l, 2, 4, 5} = J2)i3) ji}- For brevity, we simply write (J1J2J3J4) below. Because the 
contour plot at the left of Figure 10 is topologically identical to G_(D), all of these Pliicker coordinates 
A/ b correspond to dominant exponentials in the contour plot. 
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The diagram at the right shows the values of the corresponding Pliicker coordinates when we choose all 
Pj = 1 (regardless of the choice of the rrij parameters). Since only the Pliicker coordinate A 12 ,6,8(^4) = 
— 1 is negative, the singular line-solitons in the contour plot are precisely those at the boundary of the 
corresponding region; these line-solitons have types [4,6], [5,8], and [2,7]-types. 

12.3. Non-uniqueness of the evolution of the contour plots for t > 0. Consider A e Sd C 
Grk.n- If the contour plot C -00(D) is topologically identical to G-(D), then the contour plot has 
almost no dependence on the parameters mj from the parameterization of Sd- This is because the 
Pliicker coordinates corresponding to the regions of C -00(D) (representing the dominant exponentials) 
are either among the collection of minors given in Theorem 5.6 (by Remark 10.5), or determined from 
these by a "two-term" Pliicker relation. Note that the minors given in Theorem 5.6 are computed in 
terms of the parameters Pi but have no dependence on the raj's. 

Therefore it is possible to choose two different points A and A' in Sd C Grk, n whose contour plots 
for a fixed Ki < . . . K n and fixed t -C are identical (up to some exponentially small difference); we use 
the same parameters Pi but different parameters mj for defining A and A' . However, as t increases, 
those contour plots may evolve to give different patterns. 

Consider the Deodhar stratum Sd C Gr 2 , 4 , corresponding to 

w = S2S3S1S2 and v = s 2 \ls 2 . 
The Go-diagram and labeled Go-diagram are given by 
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P3 


P2 
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The matrix g is calculated as g = s 2 y 3 (p 2 )yi(pz)x 2 (m)s 2 1 , and its projection to Gr 2A is 



_ /-p 3 -m 1 0\ 
V P2 ij- 

The T-function is then given by 

TA = -(P2P3-El,2 + P3-El,4 + mE 2 .4 + P2#2,3 - E 3A ), 

where Eij := (kj — Ki) exp(#i + 0j). The contour plots of the solutions with m = and m/0 are the 
same (except for some exponentially small difference) when («0. In both cases, the plot consists of 
two line-solitons forming an X-crossing, where the parts of those solitons adjacent to the region with 
dominant exponential S3 4 (i.e. for x ^> 0) are singular, see the left of Figure 15. 

On the other hand, for t 3> 0, the contour plot with to = is topologically the same as it was for 
t <C 0, while the contour plot with in ^ has a box with dominant exponential E 2 ,4, surrounded by 
four bounded solitons (some of which are singular). See the middle and right of Figure 15. So not only 
the contour plots but also the soliton graphs are different for t >• 0! 

Note that the non-uniqueness of the evolution of the contour plot (a tropical approximation) does 
not imply the non-uniqueness of the evolution of the solution of the KP equation as t changes. If one 
makes two different choices for the TOj's, the corresponding T-functions are different, but there is only 
an exponentially small difference in the corresponding contour plots (hence the topology of the contour 
plots is identical). This is particularly interesting to compare with the totally non- negative case, where 
the soliton solution can be uniquely determined by the information in the contour plot at t -C 0. For 
more details, see the results on the inverse problem in [20]. 
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Figure 15. The non- uniqueness of the evolution of the contour plots (and soliton 
graphs). The left panel shows the contour plot at t = —20 for any value of m. The 
middle panel shows the graph at t = 20 with m = 0, and the right one shows the graph 
at t = 20 with m = 1. These contour plots were made using the choice Pi = 1 for all i, 
and (ki, . . . , K4) = (—2, —1, 0, 1.5). In all of them, the region at x 3> has a positive 
sign (A3. 4 = 1) and other regions have negative signs. This means that the solitons 
adjacent to the region for x ^> are singular. 
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